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Abstract. We present one- and two-jet inclusive cross sections for 7*7 scattering 
and virtual photoproduction in ep collisions. The hard cross sections are calculated 
in next-to-leading order QCD. Soft and collinear singularities are extracted using the 
phase-space-slicing method. The initial state singularity of the virtual photon depends 
logarithmically its' virtuality. This logarithm is large and has to be absorbed into 
the parton distribution function of the virtual photon. We define for this purpose an 
MS factorization scheme similar to the real photon case. We numerically study the 
dependence of the inclusive cross sections on the transverse energies and rapidities of 
the outgoing jets and on the photon virtuality. The ratio of the resolved to the direct 
cross section in ep collisions is compared to ZEUS data. 
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1 Introduction 

The topic of this work is the structure of the virtual photon as it can be deter- 
mined in jet production in high energy collisions. In particular, we will study 
electron proton scattering as is explored at HERA and the scattering of virtual 
on real photons as is possible at e + e~ colliders. 

In the parton model [1] a hadron is thought to consist of point-like particles 
that can be identified mainly with the valence quarks. The valence quarks are 
surrounded by a sea of virtual quarks and are bound by gluons. These parti- 
cles obey the laws of Quantum Chromodynamics (QCD) which is a non-abelian 
SU(3) gauge theory [2], Hadrons are usually probed in high energy scattering ex- 
periments. These experiments involve contributions from a wide range of scales. 
An important property of QCD is its asymptotic freedom [3], which states that 
the coupling between quarks and gluons vanishes for asymptotically small dis- 
tances. Factorization theorems allow a separation of the short and the long dis- 
tance contributions of the high energy scattering (for reviews on this subject see 
[4, 5]). This permits the application of perturbative QCD to calculate the hard 
part of the cross section. The contributions from long distances are parametrized 
by the parton distribution functions (PDF's). 

'Supported by Bundesministcrium fur Forschung und Technologie, Bonn, Germany, under 
Contract 05 7 HH 92P (0), and by EEC Program Human Capital and Mobility through Network 
Physics at High Energy Colliders under Contract CHRX-CT93-0357 (DG12 COMA). 



Inclusive Single- and Dijet Rates in NLO QCD 



2 



At the HERA collider at DESY, the scattering of leptons on protons pro- 
duces jets with large transverse energies Et- The ZEUS and HI collaborations 
have observed an important fraction of events at small virtualities P 2 ~ of the 
exchanged photon [6, 7]. The lepton is only weakly deflected in these so-called 
photoproduction events, so that it escapes unobserved in the beam direction. The 
momentum spread and the slight off-shellness of the photons that are radiated 
by the lepton, is described by the Weizsacker- Williams formula [8], where the 
photons are assumed to be real (P 2 = 0). The exchange of the other electroweak 
gauge boson, Z°, is largely suppressed for photoproduction and can be neglected. 
The transverse energy Et serves as the large scale in photoproduction, which 
allows a perturbative calculation of the hard part of the scattering. In leading 
order (LO) two different processes can be identified in the hard cross section. In 
the direct interaction, the photon couples as a point-like particle to the partons 
from the proton, leading to the Compton scattering and the photon-gluon fusion 
subprocesses. In the resolved interaction, the photon acts as a source of partons, 
which can interact with the partons from the proton. The resolved photon is 
described by the photon PDF. For quasi-real photons with virtuality P 2 ~ the 
parton content is constrained reasonably well by data from deep-inelastic 7*7 
scattering [9, 10]. Both LO processes produce two outgoing jets with large Et- 
Studies of photoproduction events with two jets in the final state at HERA have 
shown that both, the direct and the resolved processes are present for photons 
with very small virtuality P 2 ~ [6, 7]. Comparisons between theoretical pre- 
dictions for dijet photoproduction rates and the data from [6, 7] have been done 
in [11, 12, 13]. 

For the comparison between the data on jet cross sections and the theoretical 
predictions in [11, 12, 13], the hard part of the scattering has been calculated 
in next-to-leading order (NLO) QCD. In NLO, one encounters inital and final 
state singularities, due to collinear and soft radiation of partons in the initial or 
final state. There are two reasons for performing NLO calculations, which are 
far more cumbersome than the LO ones. First, one wishes to reduce the unphys- 
ical scale dependences. Second, only in NLO can one sensibly implement a jet 
algorithm, which is needed for a comparison between theory and experiment. 
However, the above discussed distinction between direct and resolved photopro- 
duction becomes ambiguous in NLO. When two-jet events are observed in an 
experiment, a disposition of energy near the beam pipe of the detector in the 
forward region of the photon can be attributed either to the photon remnant of 
a resolved photon or to a collinear final state particle from the direct interaction. 
The collinear particle in the NLO direct cross section produces a large contribu- 
tion that has to be subtracted and combined with the LO resolved term. This 
introduces a dependence of the photon PDF on the factorization scale M 1 . The 
factorization scale determines the part of the NLO direct contribution, which 
has to be absorbed into the resolved contribution. The M 7 dependences of the 
remaining NLO direct and the LO resolved contribution coming from the pho- 
ton PDF cancel to a large extent. This cancellation has been demonstrated and 
analyzed numerically in [13, 14, 15] for real photoproduction (sec [16] for related 
work) . 

Information complementary to the ep collision experiments from HERA can 
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be obtained from e + e~ colliders. Assuming the two leptons to emit quasi-real 
photons that are both described by the Weizsacker- Williams approximation, one 
effectively has 77 scattering. Both photons can be point-like or act as a source 
of partons. Three cases can be distinguished, according to the different contribu- 
tions to the cross section [17]. The interaction of a direct with a resolved photon 
is denoted as the single resolved (SR) contribution. Interactions, where both 
photons are resolved are called double-resolved (DR) contributions. These two 
cases are also encountered in ep scattering, where the one resolved photon has to 
be substituted by the proton. In addition to these possibilities, also both photons 
can interact directly in 77 scattering, which gives the direct (D) contribution. 
The region of high center of mass energies is of special interest for obtaining 
information beyond the low E T region that is determined by soft physics. This 
has been measured at LEP [18] and TRISTAN [19]. Comparison of the data in 
[19] with theoretical predictions using similar methods as those employed in [14] 
can be found in [13, 20, 21] for real photons. 

Recently, data has been presented by the ZEUS [23] and the HI [24] col- 
laborations for electron-proton collisions involving photons with small, but not- 
vanishing P 2 that allow a test of the virtual photon structure. So far, there 
has only been one measurement of the virtual photon structure function from 
the PLUTO collaboration at the PETRA e+e" collider [25]. In [26] we made a 
comparison of theoretical NLO predictions for j*p inclusive jet production with 
data from [23], by extending the methods used in [11, 12]. This extension will be 
described in detail in this work. Some theoretical studies of inclusive j*p cross 
sections in LO have been presented in [27, 28, 29]. We also include the case of 
7*7 scattering that will become important at LEP2 [30], which is an extension 
of the work from [20, 21]. 

Since the partonic subprocesses of ep and 7*7 scattering are very similar, we 
will take over the notation 'SR' and 'DR' from the 77 case to ep scattering to 
simplify the discussion. In ep scattering, the SR component denotes the contri- 
bution, where the virtual photon is directly interacting with the partons from 
the proton, whereas in the DR component the resolved photon interacts with 
the partons from the proton. 

The extension from real to virtual photoproduction is done by taking the 
Weizsacker- Williams formula to describe the momentum spread of the virtual 
photon, but keeping P 2 fixed, not integrating over the region of small P 2 and 
not assuming P 2 = 0. This is described by the unintegrated Weizsacker- Williams 
formula. In the hard process, the matrix elements for finite P 2 have to be taken. 
The matrix elements and the initial and final state singularities for the SR con- 
tribution with P 2 ^ have been calculated in [31, 32] in connection with deep- 
inelastic ep scattering (DIS) at HERA, where P 2 is large. Since we consider P 2 
to be finite, the photon initial state singularities encountered in real photopro- 
duction do not occur. Instead, when integrating over the phase-space of the final 
state particles, a logarithm of the type ln(P 2 /E^) occurs. In DIS this logarithm 
is small, since P 2 is of the order of E\, and thus has not to be considered sep- 
arately. In virtual photoproduction though, P 2 is small and the logarithm gives 
a large contribution. This large term has to be subtracted from the SR hard 
cross section, where the virtual photon is direct, and combined with the resolved 
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virtual photon from the DR contribution. This introduces a dependence of the 
virtual photon PDF on the factorization scale M 7 , just as in the case of real 
photons. The cancellation of the M 7 scale dependences of the NLO direct and 
resolved contributions must hold also for virtual photoproduction with P 2 ^ 0. 
This has been worked out in [26] and will be studied numerically in this work. 

The D contribution is needed for the direct interaction of one real and one 
virtual photon in 7*7 reactions. The initial and final state corrections for the D 
contribution are calculated here for the first time. For the real photon, the singu- 
larities are handled as discussed in [13, 21]. For the virtual photon, the procedure 
is equivalent to the one described for the SR contribution in ep scattering. 

The theoretical calculation of the resolved cross sections requires the parton 
distribution functions of real and virtual photons. Several parametrizations of 
the parton contents of the real photon are available in the literature by now 
[33, 34, 35] and seem to be consistent with dijet production data in ep scattering 
[13, 14]. For virtual photons theoretical models have been constructed that de- 
scribe the evolution with the scale Q 2 of the parton distributions and the input 
distributions at scale Qo with changing P 2 [27, 36, 37]. However, these virtual 
photon PDF's are not available in a form that parametrizes the Q 2 evolution in 
NLO. Only the LO parametrizations for the virtual photon are given in [36, 37]. 

The outline of this work is as follows. In section 2 we will discuss the general 
structure of factorization and renormalization for NLO corrections. Section 3 
contains a calculation of the LO and NLO partonic cross sections for the D, SR, 
and DR contributions with a virtual photon. We recall the virtual corrections 
to the Born cross sections and present the inital and final state singularities, 
using the phase-space-slicing method. The DR contribution has been calculated 
in [13] and will be considered only briefly. The parton distribution function of the 
virtual photon is discussed in section 4. Section 5 contains numerical results for 
inclusive single- and dijet production in ep scattering. Several numerical tests will 
be presented and the available data is compared with our theoretical predictions. 
Section 6 gives theoretical results for 7*7 collisions with the kinematics of LEP2. 
Finally we present a summary and an outlook in section 7. The appendix contains 
the analytic results for virtual, initial, and final state corrections. 

2 General Structure of the Hadronic Cross Sections 

The key to using perturbative QCD is the idea of factorization. It states that 
a cross section is a convolution of different factors that each depend only on 
physics relevant at one momentum scale. In this section we explain how factor- 
ization shows up in the hadronic cross sections we use in this work. Especially 
we will discuss the divergences appearing in a NLO calculation of the pertur- 
bative hard cross section and explain how these divergences are factorized and 
absorbed by a redefinition of the PDF's involved in each process. The general 
procedure described in this section will be applied to the specific partonic NLO 
cross sections that are calculated in section 3. 
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2.1 Factorization of Hard and Soft Regions 

The physical cross sections considered in this work have a general structure, 
where the long-distance and short-distance parts are separated. The hadronic 
cross section da H of a process is given by a convolution of the hard cross section 
do ab and the PDF's f a /A{x a ) and f b / B (x b ): 

dcr H J dx a dx b f a/A (x a )da ab f b/B (x b ) . (1) 

a. b 

In general, the PDF /^(x) of a hadron A gives the probability of finding a 
parton i (quark or gluon) with momentum fraction x within the hadron. It cannot 
be calculated perturbatively and has to be fixed by measurement. The partons 
from the PDF are thought to interact in a hard process involving a large scale 
that allows to make use of the asymptotic freedom of QCD [3], i.e. the vanishing 
of the coupling between the partons for asymptotically small distances. For a 
large scale (i, the QCD coupling constant g(fi) behaves as g(fi) ~ 1/ \n{fjb/ Aqcd) 
and a perturbative expansion of the hard cross section in the strong coupling 
constant can be applied. The hard process in (1) is described by the partonic 
cross section 

da ab = -^—\M ab \ 2 dPS^ , (2) 

ZXaX b S 

where 2x a x b s is the flux factor, |.M a {,| 2 are the partonic matrix elements and 
dPS^ represents the phase space of the n final state particles of the subprocess. 
In the final state we are interested in jets for which suitable jet definitions have 
to be defined in order to go from the partonic level to observable quantities. We 
will come back to this in section 5. The general structure of the cross sections 
discussed in this work is indicated in Fig. 1. For the case of ep scattering, A 




(P t ) 



Fig. 1. Factorization of hard and soft processes in the hadronic cross section. 

will be a lepton, that radiates a virtual photon and B will be a proton. The 
remnant A will stem from the resolved virtual photon, whereas the remnant B 
comes from the proton. In the case of 7*7 scattering, A is a virtual and B a real 
photon. They can both have a hadronic structure, leading to the remnants A 
and B. The different subprocesses, encountered in these two cases are explained 
in section 3. 
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2.2 NLO Corrections 

The matrix elements of the partonic cross section (2) can be calculated by sum- 
ming up all Feynman diagrams to a given order. We will be interested in processes 
with two initial and at least two final state particles. The calculation of the LO 
contributions is straightforward. In NLO several difficulties are encountered. We 
have to distinguish the virtual corrections to the 2-^2 partonic processes, which 
contain self energy and vertex corrections, and the real corrections, which stem 
from the radiation of an additional real parton from the 2^2 processes, leading 
to 2 — > 3 processes. Both these contributions contain characteristic divergences. 

As an example for a 2 — > 3 process the Feynman diagrams for the amplitudes 
of the Born subprocess qq — ► qq and the 0(a s ) correction containing a real gluon 
emission are drawn in Fig. 2. We consider a parton with momentum p a emitted 




Fig. 2. Partonic cross section: Born graph and real gluon emission. 

from a resolved photon. Taking k to be the momentum of the outgoing gluon, 
the Feynman diagram contains a propagator of the form 

In the limit of massless quarks, the propagator diverges in certain regions of 
phase space. The denominator 

M 2 = ( Pa -k) 2 = -2 Pa k = -2|p ||fc|(l- cos 0) , (4) 

where 9 is the angle between the gluon and the parton, vanishes if cos# = 1 
(collinear divergence) and if |fc| = (soft divergence). Both, the collinear and 
soft divergences are infra-red (IR) divergences that can be regularized in the 
dimensional regularization scheme [38, 39]. In this scheme n — 4 — 2e dimensions 
are chosen for the phase space integration, so that the singularities appear in 
poles like 1/e and 1/e 2 . After the poles have been removed the limit e — > 
is taken and the four-dimensional result is obtained. The singularity shown in 
Fig. 2 is due to the radiation of a gluon in the initial state and is thus called 
an initial state singularity. In addition to these, final state singularities occur, 
when a parton is collinear or soft in the final state. 

The virtual corrections involve loop integrals over internal momenta, that 
lead to ultra-violet (UV) and IR divergences. These divergences can be extracted, 
as the real corrections, in the dimensional regularization scheme as poles in 1/e 
and 1/e 2 . The UV divergences are removed completely by adding a counter term 
to the QCD Lagrangian, where the singularities are absorbed by a renormaliza- 
tion of the quark charge, quark field, and gluon field. 
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When the virtual and the real contributions are added, the IR divergences 
cancel partly, leaving only initial state singularities. It can be proven that in the 
remaining hard cross section the short distance finite parts and the long distance 
singular parts factorize [4, 5]. One defines a bare partonic cross section da, that 
is calculable in perturbative QCD, a renormalized finite partonic cross section 
da and transition functions so that [40] 

da ij (s)=y~] / dz a dz b r i ^ k (z a ,fi A )daki{z a z b s,HA,IJ>B)r j ^i(z b ,iJ,B) ■ (5) 
k,l 

The variables z a , z b G [0, 1] give the momentum fraction of p a ,Pb in the propa- 
gator after a parton is radiated as can be seen in Fig. 2. The singular terms are 
absorbed into the transition functions in such a way that the renormalized cross 
section is finite. This absorption depends on the scales \x A and [Ib, which are the 
factorization scales for the hadrons A and B, respectively. To obtain a hadronic 
cross section, which is free of divergent parts, one needs renormalized PDF's /, 
which are defined as 

l l 

fiA{Va,HA) = J J dxdzf jA {x)r i ^ 3 {z a ,ii A )8{rj a - x a z a ) 

o o 
l 

= f^fiA^Ar^z*,^) . (6) 

J z a \z a J 

As one sees, the factorization scale dependence of the transition functions leads 
to a scale dependence of the renormalized PDF's. The factorization of the hard 
and soft parts in the partonic cross section is pictured in Fig. 3 for the case 
of a resolved photon with the subprocess depicted in Fig. 2. The factorization 




Fig. 3. The factorization theorem for the singular part of the partonic cross 
section. 

scales [i a and \ib define what is to be understood as the hard and the soft 
part of the cross section. Referring to equation (3), M 2 gives the off-shellness of 
the propagator in the initial state of the partonic cross section. Interactions with 
M 2 < ii a are described with help of the PDF of hadron A, whereas for M 2 > \ip 
one can apply perturbative QCD and calculate the partonic cross section. 

Using the definitions of the renormalized quantities, the IR safe hadronic 
cross section reads 

da H (s) d Vadr] b f kA (r]a,VA)daki(r]aVbS,LiA,VB)fiB(Vb,VB) , (7) 

k,l J 
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where the variables rj a ,r)b G [0,1] are defined as rj a — x a z a and r\b — x^z^. 
The connection between the IR safe and the bare hadronic cross sections can 
be easily seen, by inserting the definition of the renormalized PDF's (6) into 
(7) and performing the integrations over the delta functions, making use of the 
definition (5). The factorization scale dependences of the renormalized partonic 
cross section and the PDF's cancel to a large extend. The transition functions 
that connect the renormalized partonic cross section and the PDF arc, however, 
not unique in NLO and arbitrary finite parts can be shifted from the PDF's to the 
renormalized partonic cross section. Therefore one has to define a factorization 
scheme to be used for a consistent calculation. Commonly used schemes are the 
DIS [41] and the MS [42] schemes. 

To extract the renormalized from the unrenormalized quantities, one assumes 
da, da and the transition functions to have perturbative expansions in a s [40] : 

oo 

d ^ = E(£)"^ (n) ( s ) ( g ) 

n 

da(s^) = g(g)" d(T (")( S ) (9) 

n 

oo 

r^ k (z) = Slk s(i-z) + Y / (^) n r^ k (z) (io) 

n=l 

The LO contributions are understood to be the n = 0th order contributions. For 
the DR partonic cross section it is of order 0{a 2 s ), for the SR contribution is 
of order (D(aa s ) and for the D contribution it is of order 0(a 2 ). Inserting the 
expansions (8)-(10) into (5) gives up to 0(a s ) 



2tt 



+ 



+ 



dalj'(s,fi A ,fi B ) 
E / dz 1 r i ^ k (z 1 ,f, A )da^(z 1 s) 

k 

£ J dz 2 d-a^{z 2S )r^{z 2 ^ B ) 



(11) 



Comparing the left hand and the right hand side in LO gives da^ — da^°\ The 
NLO correction is obtained by comparing the left hand and right hand side to 
order a s and rearranging the terms: 

- < - E / - E / • ( 12 ) 

k j k j 

Thus the prescription for subtracting the singular parts from the bare cross 
section is simple: the singularities are removed by a convolution of the finite 
Born cross section with the singular 0(a s ) transition functions. 

As we have seen, the PDF's acquire a dependence on the factorization scales. 
The evolution of the PDF's with the scale are predicted in perturbative QCD 
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by the DGLAP evolution equations [43, 44, 45]. This is the following set of 
integro-differential equations: 



Here, Q 2 is a general scale, and P^ (z) are the splitting functions that represent 
a process in which a parton with momentum fraction x radiates a parton with 
momentum fraction (1 — z)x and continues with momentum fraction xz. The 
splitting functions can be expanded in powers of the strong coupling constant 
with a s = g 2 /4:Tr: 



The evolution equations (13) predict the PDF's at a higher scale once they are 
fixed at some input scale Qq. 

2.3 Factorization for the Photon 

In the previous section we have described, how the absorption of the singularities 
works in the case, when A and B are hadrons or resolved photons. For the direct 
photon one additional complication has to be taken into account. 

As mentioned in the introduction, in LO the photon gives rise to direct and 
resolved contributions in the hadronic cross section. In NLO, the creation of a 
collinear qq-p&ir in the initial state leads to initial state singularities in the direct 
contribution. The only place for these singularities to be absorbed is the PDF 
of the resolved photon. This leads to a point-like term in evolution equations of 
the the photon PDF [46, 47]. A subtraction procedure for the real photon that 
is consistent with the evolution of the photon PDF has been worked out in [48] . 

For the virtual photon, one actually has no real singularity, since the virtual- 
ly P 2 regularizes the divergence. Integrating over the phase space of the gg-pair 
in the initial state leads to a logarithmic dependence on P 2 , namely \n(P 2 /Q 2 ), 
where Q 2 is the hard scale of the process. This logarithm becomes large for 
P 2 <C Q 2 and is absorbed into the PDF in much the same way as described in 
[48]. This leads to an inhomogcneous term in the PDF of the virtual photon, 
which differs somewhat from the point-like term in the case of the real photon. 
This will be described in more detail in section 4, where the construction of the 
virtual photon PDF is explained. 

3 Partonic Cross Sections 

In this more technical section, we proceed with a computation of the perturba- 
tively calculable partonic cross sections. The partonic cross sections contributing 
to ep and 7*7 scattering are very similar and will therefore be treated together 
in this section. Both hadronic cross sections contain single resolved (SR) con- 
tributions, in which the virtual photon couples directly to the subprocess and 



1 




X 



(13) 




(14) 
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double resolved (DR) contributions, in which the virtual photon is resolved. In ep 
scattering the virtual photon and its partonic content interact with the partons 
of the proton, whereas in 7*7 scattering they interact with the parton content 
of the real photon. In addition to the SR and DR contribution, in 7*7 scattering 
one encounters the direct (D) interaction of both photons. 

After an introduction to the notation of the various relevant subprocesses, 
we give the formula? for the Born and virtual contributions. Then we explain the 
phase-space slicing method as a tool to separate singular regions of phase space in 
the partonic cross section, so that we can calculate the singular parts of the real 
final and initial state corrections. The results from [31] for the SR contributions 
are recalled for completeness and consistency. The DR contributions can be 
found in [13] and will be considered only briefly. 

3.1 Notation 

For the calculation of hadronic cross sections in sections 5 and 6, we will have 
to calculate the matrix elements for the various partonic cross sections in LO 
and in NLO. In LO the D contribution is of 0(a 2 ), the SR contribution is 
of 0(aa s ) and the DR contribution is of 0(a 2 s ). In NLO the D, SR and DR 
contributions are of one order higher in a s . Since the partonic cross sections 
have to be convoluted with the PDF's of the photon, which are of order a/a s in 
the high energy limit [46], the different contributions will turn out to be of the 
same order. The matrix elements \M\ 2 from equation (2) are obtained by taking 
the trace of the hadron tensor that corresponds to each subprocess. We define 
H = —g IMV H liV . The Born contributions are labeled H Bl the virtual corrections 
are Hy, and the real corrections are H B . 

First, we have collected the definition of the LO Born matrix elements in 
Tab. 1. The matrix elements for incoming anti-quarks are the same as those for 
quarks and only give a factor 2 in the sum over all contributions in the hadronic 
cross section. This holds also for the contributions in the other tables. Of special 
interest in this work will be the NLO corrections to the Born matrix elements for 
processes involving a virtual photon. These are the D and SR contributions, that 
are collected in Tab. 2. The SR contributions with one virtual photon have been 



Table 1. Definition of the LO matrix elements. 



D and SR Processes 


DR Contributions 


Bi = i?s(7*7 — > qq) 


B i = H B {qq l ^qq') 


B 9 = H B (qq — > gg) 


B 2 = H B (~f*q -> qg) 


B 5 = H B (qq' - qq') 


Bio = H B (qg qg) 


B 3 = H B <a*g -> qq) 


Be = H B (qq ->■ qq) 


Bn = H B (gg -> qq) 




B r = H B (qq^q'q') 


B12 = H B (gg -> gg) 




B 8 = H B (qq -> qq) 
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Table 2. Definition of the NLO matrix elements for the D and SR processes. 



Virtual Corrections 2^2 


Real Corrections 2^3 


= JT V (7*7 - ««) 


#i = ^(7*7 -» qqg) 


i?4 = H R (^*q -> qqq) 


V 2 = H v (j*q -> qg) 


H 2 = H R (-f*q -> 93.9) 


H 5 = H R ("f*g -> qqg) 


V 3 = H v (Y~g -> 99~) 


#3 = fffl(7*9 - MY) 





Table 3. Definition of the NLO matrix elements for the DR contribution. 



Virtual Corrections 2^2 


Real Corrections 2^3 


V 4 = H v (qq' 


#6 = H R (qq' -> ggr'p) 


#7 = H R (qq' -> g<fg) 


V 5 = H v {qq> - qq') 


#8 = #fl(99 995) 


#9 = iifl(99 -» 9'9'ff) 


V 6 = H v (qq -> qq) 


#10 = #fl(99 -» 99.9) 


#11 = #r(9.9 99' 9') 


V 7 = H v {qq^q'q') 


-ffi2 = H R (qg — > ggg) 


#13 = H R (qq -> g.gg) 


V 8 = H v (qq -> g<?) 


i?i4 = #k(9.9 -»■ 953) 


#15 = #k(3.9 99ff) 


V 9 = H v (qq-> gg) 


#16 = #it(53 -»■ 553) 




V w = H v (qg -> gg) 






Vii = i?y(gg -*■ 99~) 






V12 = #v(33 -> 55) 







studied by several authors [31, 32], the D and SR contributions for real photons 
can be found in [13], whereas the D contributions with one virtual photon have 
only been studied in this work, yet. 

Because of their importance, we show in Fig. 4 the classes of matrix elements 
as collected in Tab. 2 explicitly. We also show the definition of the momenta 
which will be used throughout this work for all three, i.e. D, SR and DR, con- 
tributions. In Tab. 3 the matrix elements of the NLO processes for the DR 
case are collected. These matrix elements have been calculated in [49] and the 
integrations over the singular regions of phase space where performed in [13]. 

3.2 The Two-Body Processes 

For the 2^2 processes we use the Mandelstam variables 



S = (Pa + Pbf = (Pi +P2) 2 
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Fig. 4. Notation of the different processes involving one virtual photon. 



t = (Pa-Pi) 2 = {Pb-P2? , (15) 

U = (Pa -P2) 2 = (Pb -Plf ■ 

Note, that in the D and SR case, p a = q with P 2 = —q 2 . In the D case pt is a 
real photon, in the SR and DR cases pi, is a massless parton. The partonic cross 
section is given by the flux factor, the two-particle phase space and the matrix 
elements: 

da(ij -jets) = li/fe-^jcts)dPS (2) . (16) 
The two-particle phase space is given by 

where z = {pbPi) / (PaPb) ■ Expressed by the Mandelstam variables, the phase 
space reads 

d PS( 2 ) = _l_filY( s + P 2 )- 1+2e - , (18) 
r(l-e) \stuj K ' 8tt ' v ' 

if the particle p a = q has mass — P 2 . This is valid for the D and SR case. In the 
DR case all partons are massless, so we substitute P 2 = and the phase space 
reduces to 

^ = nhrM'£. ■ <19) 

The Born matrix elements for the D and SR case read, using the notation of 
Tab. 1, 

B x = (16^ V) (Qt8N c ) T 7 (s,i,u) , (20) 
B 2 = -(16ir 2 aa s ) (Q 2 2C F ) T 7 (u,i,s) , (21) 
B 3 = (l6TT 2 aa s ) (Q 2 ) T 7 (s,t,u) , (22) 
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where the definition of T 7 (s, t, u) can be found in the appendix, section 8.1. The 
Born matrix elements for the DR case are given by 





= (167r 2 a 2 ) -r±T Tds.t.u) , 




(23) 


Bk 


= (167T 2 a 2 ) -rkrr Ti (u t s) 




(24) 


B 6 


= (l6TT 2 a 2 s ) \ [Ti(8,t,u)+Ti{8, 


u,t)+T 2 («,t,u)] , 


(25) 


B 7 


= (16?r 2 a 2 ) Ti («,*,*) , 




(26) 


B 8 


= (167r 2 a 2 s ) j^t [T 1 (u,t,s)+T 1 (u,s 


,t)+T 2 (u,*,«)] , 


(27) 


B 9 


= (l6TT 2 a 2 s ) T 3 (s,t,u) , 




(28) 


Bio 


= -(167r 2 a 2 ) T 3 (i, S , U ) , 




(29) 


B u 


= ( 167r2a s) 16n\ Cf T 3(s,t,u) , 




(30) 


B\2 


- (167r 2 a 2 ) 2i 6 jv^c F r 4(s,t,u) . 




(31) 



The factors 1/2 in some of the expressions are symmetry factors for two identical 
particles in the final state. The definitions of the matrix elements Ti, . . . , T4 can 
again be found in the appendix 8.1. 

The virtual corrections for the SR case are calculated by multiplying the one- 
loop diagrams for 7*7 — > qq, -f*q — ► gq and 7*5 — ► qq with the corresponding 
Born diagrams, which leads to an extra factor a s in the matrix elements. The 
virtual corrections V2 and V3 are well known for quite some years from the 
processes e + e~ — > qqg [53, 54]. They are achieved by crossing from the known 
matrix elements [31]. The contribution V\ for the D case can be inferred from the 
SR case by considering the contribution H v (j*g — > qq)- Only the parts which 
have no gluon self-coupling are taken into account and the color factors have 
to be adjusted appropriately. To compare the singular structure of the virtual 
corrections with those from the real corrections, we write down the final result 
for the D and SR case in the form 



Vi 



V 2 = 



IQtt z 



~2~r 

™,2 



r(l-e) 



r(i 

2 



l 6 ^_£ 2( l_ e) 



1 (\ 



2e) 

r(i 



<0 



r(l-2e) 



V 3 = 16tt 2 



- ! - Y, - 11 A 

47T/X 2 



2?r 



C F T 7 (s,i,u) 
r(l-e) 
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r(l - 2e) 
T 7 (s, u, t) 



CfE 3 



C F E\ — -NcCfE2 



-N C E A 



(32) 



(33) 



(34) 



Terms of order 0(e) have been neglected. The expressions , . . . , E4 are given 
in the appendix 8.2. In the DR case, the Born processes B4, . . . , B 12 have to be 
multiplied by the corresponding one-loop processes. The results can be found in 
[13]. They are given by 14, ... , V 12 in [22]. 
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3.3 Phase-Space-Slicing Method 

As discussed in the section 2, the partonic 2^3 corrections are singular in cer- 
tain regions of phase space. One possibility is to integrate the hard cross section 
over the complete phase space in n dimensions and to obtain analytical results 
for the integrated matrix elements. We choose a somewhat different method in 
this work. We separate the singular regions of phase space from the finite regions 
by inserting an invariant mass cutt-off y c into the integration, symbolically 

1 Vc l 

J dPS (3) |M 2 ^ 3 | 2 = J dPS (3) |M 2 ^ 3 | 2 + J dPS (3) |A4 2 ^ 3 | 2 . (35) 

y c 

The first integral on the right hand side of this equation contains the singular 
region of phase space and is integrated analytically in n = 4 — 2e space time 
dimensions. If the cut-off y c is chosen appropriately small, this singular phase 
space separates into a 2 — > 2 phase space dPS^ that will be kept and a remaining 
part dPS^* 1 that is integrated out together with the matrix elements. For small 
y c various approximations can be applied to the matrix elements, so that the 
integration of the matrix elements over e/PS^* 1 is simplified considerably as will 
become evident in sections 3.4-3.6. The non-singular second term on the right 
hand side of (35) is integrated numerically, opening the possibility to adopt 
a wide range of jet definitions and experimental cuts. This flexibility allows a 
detailed comparison between theory and experiment. Of course, the dependences 
of the first and the second contribution on the parameter y c should compensate, 
leaving a result independent of y c , since the cut-off has no physical significance. 
This also provides a strong test of the results, which will be described in section 
5. The method described here is referred to as the phase-space-slicing (PSS) 
method [50]. 

An important step in the application of the PSS method is the separation of 
the two-body phase space from the singular part, c£PS^ — > dPS^ciPS^. This 
separation is different for the three cases encountered in this work, which are 
singularities in the final state, in the initial state for massless particles and in 
the initial state for a massive virtual photon. In the next two subsections we 
will provide the formulae that serve as a basis for the calculation of the singular 
parts of the partonic cross sections in sections 3.4-3.6. 



3.3.1 Singularities in the Final and in the Massless Initial State We 

consider the splitting for a phase space containing one massive particle p a with 
mass P 2 = —p 2 ,. In general, the 2^3 phase space in n dimensions is given by 
[51] 



dFS = 2^(2^-1 2^(2.)^ 2E 3 (2n)^ ^ * + 

(36) 



It is useful to introduce the following irreducible set of invariants: 



,s = ZpaPb - P 2 , h = -2p b pi, 
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si=2p lP2 , t 2 = -2 PaP3 -P 2 , (37) 

s 2 = 2p 2 p 3 . 

These five invariants are pictured in Fig. 5 on the left. The separation pictured 
on the right of Fig. 5 is achieved by inserting [52] 

into (36), where the definition s 2 = p 23 = (p 2 +P3) 2 is used and E 23 is the energy 
of this intermediate particle. One obtains 

( d n ^ 1 r> d n ^ 1 n 1 

X \ 2E 2{ 2^ 2E , {2 :;-^ -P,-Ps)(2,r} . (39) 

To perform the integration over the delta functions in (39) it is useful to define 
the kinematical variables in the cm. system of the outgoing partons p 2 and ^3. 
The angles of the partons p\ and p 2 ,p 3 with respect to the parton pb are shown 
in Fig. 6. To parametrize the angles, the variables b = 5(1 — cos#) and z\ = 
(PbPi) I (PaPb) are used. After integrating over the delta functions and expressing 
the variables by z\, b and the irreducible invariants introduced in (37), the three 
particle phase space in n = 4 — 2e dimensions reads 



-2e 



dFb " 2ir M^N b [b{l b)l N^r(2-2e) 



16tt 2 



s 2 zi(s (l - z\) - s 2 ) 



(40) 



with the normalization constants 
1 

N b = J db[b(l-b)}- £ = r(2 '_ £ and N4 = Jdo*m 




9 r 2 (i-e) 

(41) 

For singularities in the final state the invariants s 2 and b will vanish, whereas 
for the singularities in the initial state of the massless particle p b the invariant 
ti and therefore z\ will vanish. 




Fig. 5. Separation of the three particle phase space. 
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Fig. 6. The three particle final state in the cm. system of the particles p 2 and 

P3- 



3.3.2 Singularities in the Initial State of a Virtual Photon For this 
case, the mass P 2 = —p 2 of the photon with momentum p a serves as a regulator 
for the integration. Therefore, the phase space (36) can be calculated in n = 4 
dimensions. We introduce the five invariants 



s = 2p a p b - 
si = 2pip 2 , 
s 2 = 2p 2 p 3 , 



pz 



tl = 

*2 = 



~2pbP3, 
-2p a P\ " 



(42) 



and separate the phase space analogously to the case discussed above by inserting 
(38) into (36). Again we move to the cm. frame of the particles p 2 and p 3 and 
define the variable b = ^(1 — cos#), but now 6 is defined as the angle between p a 
and p 3 . The singularities occur for t 2 — > 0, which we parametrize by the variable 
z 2 = {PaPs) I {PaPb)- It is obvious from the definition of t 2 that the variable z 2 
vanishes only for P 2 — > 0. After integrating over the delta functions, the result 
is simply 

7T 2tt 8tt 8tt y ' 



3.4 Final State Singularities 

The singularities in the final state appear when the invariant s 2 in equation (37) 
becomes on- mass-shell. We define the variable r = s 2 /sq and consider the limit 
r — > 0. We start by considering the D and SR cases, for which we define the 
two-body Mandelstam variables as 



s = (p b + q) 2 = 2 Pb q-P 2 , 
t = {p b - Pif = -2p b pi . 

U = (Pb-P23f = -2p b p 23 



(44) 
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In the limit r — > the five invariants (37) reduce to these variables, with Sq — ► s 
and t\ — > t. The definitions of s, t and u are only unique in the limit r — ► 0. In the 
limit s 2 — > the phase space (40) separates according to dPS^ 3 ' = dPS^ctPS^ 
with 



where 



G F (r) 



1 



1 + O(r) . (46) 



The two-body phase space dPS*- 2 ** is given by equation (18). The limits of in- 
tegration in dPS (r) are given by r £ [0, -t/(s + P 2 )],b £ [0, 1] and G [0, 7r]. 
The invariant S2 is integrated up to S2 < y c so, which restricts the range of r to 
< r < min[-i/(s + P 2 ), y c ] = y F . 

We have now achieved the separation of the phase space and have to inte- 
grate the matrix elements over the region dPS (r) . Therefore the 2 -> 3 matrix 
elements H\, . . . H 5 are expressed by the variables s, t, u, r, b and in the limit 
r — > 0, which leads to approximated matrix elements with final state singulari- 
ties H F i, . . . , H F5 [31]. A difficulty arises for those squared matrix elements that 
contain real gluons. In that case more than one invariant can vanish in a propa- 
gator, so that the different classes of singularities, such as initial and final state 
singularities, are not properly separated. The separation is achieved by partially 
fractioning the matrix elements. 

After the matrix elements have been partially fractioned and approximated 
in the limit r —> 0, the integration over the singular region of phase space yields 



/ 



dPS^H Fi = Sn(^y ^-^(l-e)F i + 0(e) . (47) 



The final results for Fi, . . . , F$ are listed in appendix 8.3. They contain the IR 
collinear and soft singularities that cancel against those of the virtual corrections. 
It is essential that the singular terms are proportional to the LO matrix elements 
T 1 and that the variables s, t and u defined in (44) correspond to the two-body 
variables in the above discussed limit. The contributions F 2 , . . . , F 5 have been 
calculated in [31]. The final state correction i*\ of the D contribution can be 
derived from F 5 by keeping only the abelian part and adjusting the color factor. 
One can compare the singularities for the case of direct real photons as stated in 
[22] with the ones given in appendix 8.3. The real photon contributions F\ , . . . , F$ 
in [22] follow from the contributions F\, . . . , F$ in this work by taking the limit 
P 2 — > 0, so that t + u = —s. 

Turning to the resolved photon case, the phase space is obtained from the 
formula (45) by substituting q — > p a with p 2 = 0, so that s = 2p a pt,. The 
two-body variables in the resolved case are given by 

s = 2p a p b , t=-2p b p 1 , u = -2p b p 23 . (48) 

Note, that the two-body phase space for the DR case is given by equation 
(19) of section 3.2. Expressing the matrix elements for the resolved processes 
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Hq, . . . , Hiq as classified in Tab. 3 through the variables s, t, u, r, b and (f> yields 
the matrix elements containing final state singularities Hfq, . . . , Hpi§ [13, 22]. 
The integral over the singular phase space is the same as in equation (47). The 
final state corrections for the resolved case F 6 , . . . , Fi 6 will not be stated here 
again since they can be found in [22] . 

3.5 Initial State Singularities for Massless Particles 

We turn to the discussion of the initial state singularities for particles with zero 
mass. This includes the photon initial state singularity for the real photon for 
the D case, the parton initial state singularities in the SR case and the parton 
initial state singularities in the DR case. We start with the SR case, from which 
the others can be inferred. 



3.5.1 Parton Initial State Singularities in the SR Case In the SR case 
parton initial state singularities arise for the incoming particle p b - The incoming 
particle p a is the virtual photon with p a = q and P 2 = —q 2 . Using the notation 
of section 3.3.1, the singularities appear when the invariant t\ becomes on-mass- 
shell, i.e. for z\ — > 0. The invariant does not vanish in the case of initial state 
singularities but rather defines the partonic cm. energy of the corresponding 
two-body process. We define the new variable 



Zb = 



P2P3 

qpb 



S2 



so 



P 2 



e [%,i] 



(49) 



that gives the fraction of the momentum p b that participates in the subprocess 
after a particle has been radiated in the initial state. The variable r\b is given by 
i]b — XbZb- We define the two-body variables as 



S = (P2 +P3) 2 = "2P2P3 , 
t = (z b p b - p 2 ) 2 = -2z b PbP2 
u = (z b Pb - P3) 2 = -2z b p b P3 



(50) 



In the limit z — > the variable S2 reduces to s. In the same limit the phase space 
(40) separates according to dPS (3) = rfPS (2) dPS (r) , where 



dPS (r) = 



1 



T(l - e) Nff, 

_ e dz b ( 1 
dz\z x 



sin 



rG/(*i) 



Zb 



Zb 



Zb 



(-)'■ 

\ s J 16tt 2 



P 2 (l-z b ) 



with 



G!( Zl ) 



1-2!- 



z b P 2 



= l + 0{ Zl ) 



(51) 



(52) 



- s(l - z b ) - z h P 2 _ 

The two-body phase space is given by equation (18). The integration over dPS^ 
with Z\ G [0, —u/ (s + P 2 )], z b G [rj b , 1] and (f> G [0, n] is restricted to the singular 
region of Z\ in the range < Z\ < min{— u/ (s + P 2 ),y c } = y T . 
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Expressing the matrix elements for the direct photon case, listed in Tab. 2, 
with the variables s, t, u, z\, z b , b and <fi and taking the limit z\ — > 0, one obtains 
the matrix elements Hj 2 , . . . , Hj 5 [31] that contain initial state singularities on 
the parton side B. These are integrated according to 

J^ {r) Hh = }^{^yf^(l-e)l! + 0(e) . (53) 

Vb 

where results for l\, . . . , 7| are written down in appendix 8.4. They can also be 
found in [31]. Apart from the two-body variables s,t, u and the cut-off yi, they 
still depend on the integration variable z b . The results for the l\ contain IR 
singularities proportional to 1/e 2 that cancel against the corresponding singu- 
larities in the virtual corrections. The remaining singular parts are proportional 
to 1/e and to the Altarelli-Parisi kernels in four dimensions. These are removed 
by a redefinition of the PDF's that are the source of particle p b , which can be 
a hadron or a resolved photon. The redefinition, as explained in section 2.2, is 
achieved by equation (6) for the PDF's, introducing the factorization scale (M b ) 
dependence through the transition functions : 

l 

M 2 ) (^(1 - *>) + ^rV(z b , M 2 )) f jB (|) . (54) 



Vb 



Here fjB (j%) * s the PDF of hadron B in LO before absorption of the collinear 
singularities. The NLO transition functions are given by 



r^( Zb ,Mi) = - 1 -P i ^ {Zh) L^-± ) - r, ,,,,,, 



with dj = in the MS scheme. The renormalized partonic cross section da(j*i — ► 
jets) for parton initial state singularities is calculated from the unrenormalized 
cross section da by 

da(ri - jets) = dvtfi - jets)-^ £ J dz b r^{z b , M 2 )da B ( 7 *j - jets) . 

(56) 

The da B, s denote the Born level partonic cross sections that can be found in 
section 3.2. The factor Att^/M^ in equation (55) is combined with the factor 
47r/x 2 /s in equation (53) and leads to a factorization scale dependent term of the 
form 



1 , N 

— Pi^j(z b ) 
e 



V s J \M 2 



M, 



-P^(z b )ln[ -±) . (57) 



In this way, the subtracted partonic cross section will depend on the scale M 2 , 
as does the PDF of the hadron B, fi B . 
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3.5.2 Real Photon Initial State Singularities in the D Case In the D 

case the direct real photon can split into a qq pair that gives rise to a collinear 
singularity if the partons are emitted parallel. The calculation proceeds along 
the lines that have been described for the parton initial state singularities in 
the previous section 3.5.1. The two-body variables are defined as in equations 
(50). The phase space separation as well as the formula (53) for the integration 
over the singular regions remains unchanged. For the singular matrix element 
H\ 1 = if) (7*7 — > qqg) we have 

Vb 

where l\ is stated in the appendix 8.4. As remarked for the virtual and final 
state corrections, the initial state correction for the D case can be infered from 
the SR case with an incoming gluon instead of the real photon by dropping the 
non-abelian terms and adjusting the color factor. The singularity appearing in 
l\ is proportional to the splitting function P q ^ 1 (zi ) ) given in appendix 8.1. This 
function appears in the evolution equation of the PDF of the real photon as an 
inhomogencous (so-called point-like) term, as will be explained in more detail in 
section 4. Therefore, the photon initial state singularities can be absorbed into 
the real photon PDF, according to the procedure given in [48]. We define the 
renormalized PDF f qe of a quark q in the electron as 

1 

f qe (r, b ,M*) = f ^ (6 qj 6(l z b ) + grW 7 (z b ,M b 2 )) f, e {^j . (59) 

Vb 

The NLO transition functions are given by 

rW 7 (* 6 ,M 6 2 ) = -i Pq ^ Zb )I^± {^fj + c ^ {zb) (60) 

with C ql = in the MS scheme. In the discussed order, P g ^{z) = 0. The 
partonic cross section ^(7(7*7 — > jets) for the photon initial state singularity is 
calculated from the unrenormalized cross section da by 

*r(7*7 -» jets) = da( 7 * 7 -> jets) - ^ J dz b r^(z b , M^)da B { 1 *q -» jets) . 

(61) 

The cross section da B contains the LO virtual photon-parton scattering matrix 
clement B2 given in section 3.2. The dependence of the real photon PDF on the 
factorization scale enters in the same way as discussed for the SR case, section 
3.4.1. 

3.5.3 Parton Initial State Singularities in the DR Case The calculation 
of the initial state singularities for the DR case is very similar to the calcula- 
tions shown for the SR case, only now both incoming partons are massless. The 
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calculations for the partons p a and pb yield identical results and thus we have 
to consider these singularities only once. The singularities occur in the region 
zi — > 0. The formula for the phase space separation is obtained from equation 
(51) by substituting q p a with p 2 = P 2 = 0, so that s — 2zbp a Pb- With these 
substitutions, the phase space (51) reduces to 



dPS (r) 
with 



1 d 4> --2e,V 47r \ s n i u - e dz b fl-Z b 

— — sin d> \ — ^ o W (zi )dziz, 

r(l-e)N^ \sj I6ir 2 niJ 11 z b \ z b 



Gi(Z!) = 



1 - 



Zl 



1 - z b 



= l + 0( Zl ) 



(62) 
(63) 



The two-body phase space in the case of initial state singularities is given by 
(19). Expressing the resolved matrix elements H$, . . . , H\§ with the variables 
s,t,u, z\, z b ,b and <ft in the limit z± — > leads to the resolved matrix elements 
containing initial state singularities Hj 6 , . . . , Hj le [13, 22]. These are integrated 
similar to equation (53), leading to the results Jq, . . . , Jiq. These are not stated 
here, they can be found in [13, 22]. The initial state singularities on the proton 
side Jq, . . . , Ji6 are given by Rq, . . .i?" 6 in [22]. The cancellation of the poles 
from the real and virtual corrections proceeds as in the SR case discussed above. 
The remaining poles in 1/e are proportional to the Altarclli-Parisi kernels and 
are absorbed into the PDF's of the hadrons A and B that emit the particles p a 
and pb- 



3.6 Initial State Singularities for the Virtual Photon 

The initial state singularities described in the previous section were extracted 
and handled in the dimensional regularization scheme, i.e. in d = 4 — 2e dimen- 
sions. In the case of the real photon this is necessary, because the real photon is 
massless. The singular terms are proportional to a simple pole in e multiplied by 
the splitting function P g ^ 7 . These initial state singularities are absorbed into 
the PDF of the real photon. The NLO correction for the direct virtual photon 
becomes singular only in the limit P 2 — > 0. After integrating the phase space up 
to the invariant cut-off y c , the logarithm ln(P 2 /s) will occur that becomes large 
in the limit of small P 2 . The logarithm has to be absorbed into the PDF of the 
virtual photon, instead of the 1/e poles in the real photon case. 

To show the subtraction of the logarithm explicitly, we start by defining the 
two-body variables for the virtual photon initial state singularities. They are 
given by 

S = (P2 +P3) 2 = 2p 2 P3 , 

* = (Pb-P2? = -2p bP 2 , (64) 
u = (pb - P3) 2 = -%PbP3 ■ 

We define the variable z a as 



(65) 
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with rj a = x a z a . It gives the momentum fraction of the three-body cm. energy 
that participates in the two-body process. The definition of the three-body vari- 
ables is given in section 3.3.2, equation (42). As mentioned above, the mass P 2 
regularizes the initial state singularities of the virtual photon. The singular terms 
appear in the case t 2 — > in (42), which corresponds to z 2 — > for P 2 — ► 0. For 
z 2 — > the phase space (43) separates according to dPS^ 3 - 1 = dPS^dPS^, with 

dPZ™ = ******* and dP S W = -L * . (66) 

167H 7T Z a 87T S + P Z 

The limits of integration are given by z 2 G [0, — t/(s + P 2 )], z a G [VaA] an d 
4> G [0, 7r]. Since the integration of z 2 is restricted to the singular region we 
define the integration range for z 2 by < z 2 < min{— tj (s + P 2 ),y c } = yj. 

Expressing the matrix elements for the direct photon case, listed in Tab. 2, 
with the help of the variables s, t, u, z 2 , Zb, b and 4> and taking the limit z 2 — ► 0, we 
obtain matrix elements Hf 1 , . . . , H f 5 that contain the initial state singularities 
on the virtual photon side. These are integrated according to 

1 

J dPS {r) H? t - J ^8ttI? . (67) 

The results for J", . . . , I£ are collected in appendix 8.5. They contain the sin- 
gularities for the D and SR contributions. Apart from s,t, u and yj, the results 
also depend on the integration variable z a . All five expressions if , . . . , I§ contain 
the term 

M<p3 > = ^-<=«> k ( 1 + ^) < 68 » 

which is large for P 2 -C s and singular for P 2 = as expected. The large 
contribution has to be subtracted and absorbed into the PDF of the virtual 
photon. Here, we have the same freedom as in the case of the real photon, as has 
been described above for the D case. Finite parts can be shifted from the PDF 
to the direct cross section and vice versa. 

However, the virtual photon PDF's used later on in this work are constructed 
in a scheme similar to the MS scheme for real photons and we have to use the 
same scheme to obtain consistent results. Therefore, we subtract those terms 
that will yield the MS scheme of the real photon in the limit P 2 — > 0. In order 
to make the comparison with the case of the real photon possible, we state here 
the singular parts of the expressions If, . . .if, that appear for the real photon 
in d = 4 — 2e dimensions. They are given by [14] 

M =-\>k p ^ )+ ^ p ^H^ yj ) +1 2 • (69) 

The characteristic singularity proportional to 1/e is subtracted by absorbing the 
transition function 

r^(z a ,M 2 ) = -l P9 ^ M ^Z± [^fj (70) 
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into the PDF of the real photon. This subtraction produces a factorization scale 
dependence of the photon PDF and gives the finite contributions to the cross 
section. The expression remaining after the absorption is, in the MS scheme: 



In order to obtain the same finite terms in M(P 2 ) from equation (68) in the 
limit P 2 — > for the virtual photon case, we absorb the transition function 



(M 2 \ 
P*{1- Za ) ) P ^ Za) ~ Nc (72) 

into the PDF of the virtual photon. This leaves the finite term 

= -sfe f "W ( p? + - «.) ) + \ ■ < 73) 

that reduces to the expression Mjjg in (71) for real photons in the limit P 2 — > 
0. We therefore call this form of factorization the MS factorization for virtual 
photons. 



4 Parton Distribution Function of the Photon 

As mentioned in the introduction and worked out in section 3, the photon pro- 
duces a qq-pair in the initial state of the NLO direct contribution that leads to 
a large logarithm for virtual photons and a singularity for real photons. These 
terms have to be absorbed into the PDF of the virtual and real resolved photons, 
respectively, leading to a point-like term in the evolution equations of the photon 
PDF's. In this section, we wish to introduce the PDF of the real and the virtual 
photon. 

After a general discussion of the origin of the photon structure, we define the 
structure function and the PDF of the real photon. The evolution equations will 
be explained and the differences to the proton will be pointed out. A discussion 
of the formalism for the virtual photon PDF resembling the formalism for the 
real photon follows. Finally, we compare two parametrizations of the virtual 
photon PDF which we use in our computations. 

4.1 Origin of the Photon Structure 

The photon is the elementary gauge boson of QED. However one knows from 
soft low energy reactions, that a photon can behave like a hadron. If the time 
of the interaction between the proton and the photon is much smaller than the 
fluctuation time tf of the gg-pair, the pair will interact with the proton rather 
than with the photon itself and will give rise to a hadronic structure of the 
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photon. The fluctuation time for high energy photons with virtuality P 2 can be 
estimated from the uncertainty principle by [55] 

" = ^5 • (74 » 

where q is the energy of the photon and m q q is the mass of the pair. As P 2 
increases, tf becomes smaller, giving back the photon its structureless character. 
Thus, one can identify the direct photon, which interacts directly as a structure- 
less object, and the resolved photon, which has a hadronic structure. 

At this point it is important to further distinguish the possible configurations 
of the gg-pair which yield different contributions to the resolved photon. The 
photon can create a large size, asymmetric configuration with small transverse 
momentum kr that gives rise to soft non-perturbative effects and a small size, 
symmetric configuration with large kx that yields hard, perturbative interactions 
[56] . The soft part will behave more like a hadron and it will therefore be called 
the hadronic part of the resolved photon, whereas the hard part behaves more 
like a point-like photon and will therefore be called the point-like part. A possible 
physics interpretation of the soft part of the resolved photon is the fluctuation 
of the gg-pair into a vector-meson with m q q ~ my, which is described by the 
vector- meson dominance (VMD) model [57]. The coupling of the photon to the 
vector-meson 4£ has been predicted by the VMD model, giving 

Jv 

f 2 f 2 f 2 

i_p_ . Ju . £v_ —9 : i : 7 . (75) 

47T 47T 47T 

These ratios have been confirmed by measurements of the reaction e+e~ — ► 
hadrons. 

We have introduced the distinction between direct and resolved photons, 
but this distinction is unambiguous only in LO. In NLO the direct and resolved 
parts of the photon become intermixed through the point-like part of the resolved 
photon. A possibility to distinguish between the direct and the resolved photon 
interaction has been suggested by Levy [56]. Consider photon-gluon fusion as 
shown in Fig. 7. The diagram on the left is usually denoted as a direct process 




Fig. 7. Direct and resolved photon processes for photon-gluon fusion. 

in LO. The diagram on the right describes the fluctuation of the photon into a 
qq pair with a given kr followed by the interaction of one of the partons with a 
gluon, which produces a final state with some pr- If &t -C Pt, the process can 
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be called a resolved interaction. For kr 3> Pt the pr is too small for the final 
state partons to form two separate jets, so the diagram looks like the diagram 
on the left und thus can be considered as a direct interaction. For large P 2 it is 
more likely that kr ~> Pt so that the direct component is dominant. For low P 2 
the resolved part will be more dominant. 



4.2 Parton Distribution Functions of the Real Photon 

The structure of the real photon has been analyzed in the process — > eX, 
depicted in Fig. 8. We define x = Q 2 /(Q 2 + W 2 ) with Q 2 = -q 2 and P 2 = -p 2 , 
where W is the cm. energy of the 7*7 system and y = (pq)/(pk). Denoting the 
longitudinal polarization state of the photon as I and the transversal one as t, 
we can define photon structure functions 

= & ? h"« • (76) 
= ^ {n ' +n ' } ■ (77> 

Using these definitions, the cross section for e-f scattering can be written as 
da(ej — > eX) Anas 



dxdy Q 4 



[(l-y)F2 + xy 2 F?] . (78) 



This is in complete analogy to the DIS ep reaction. The difference to the case of 
ep scattering lies in the fact that the photon structure function can be calculated 
perturbatively in the limit of large Q 2 . This is not possible for the proton struc- 
ture function. The photon structure function is computable in the quark-parton 
model (QPM) from the box diagram 7*7 — > qq and gives in LO in the limit 
to 2 ,. <C Q 2 1 where m qi are the quark masses, [59, 60, 61] 

Fl* 1 = J2xqHx,Q 2 ) (79) 

i=l 

with 



L u ■ (J' 1 :<f - \ [x 2 + (1 - x) 2 ] In Q2{ \ X) + 8s(l - x) - l) . (80) 

Z7T [ m qi X ) 




Fig. 8. Single-tag DIS ej experiment. 
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The function qj can be interpreted as the PDF of the quark in the photon, in 
analogy to the proton case. 

The QPM result can be modified substantially by QCD effects, such as mul- 
tiple gluon radiation ladder diagrams. In the limit of large Q 2 these kind of 
corrections where shown to be exactly calculable in LO [46] and NLO [47] . The 
result is of the form 



a(x) , , . 
<b(x) 



a s (Q 2 ) 



(81) 



where a(x) is the LO and b(x) the NLO result. Unfortunately F^ asymp be- 
comes negative for small values of x, which cannot be true, since the photon 
structure function is measurable. The problem cannot be cured by adding the 
VMD contributions that have been mentioned in the previous section, since this 
non-perturbative contribution is expected to be well-behaved. Therefore it is not 
possible to compute the photon structure function by perturbation theory alone. 

To handle the problems of the photon structure function, Gliick and Reya [62] 
have suggested to formally add all contributions to the photon structure, namely 
the QPM, their QCD corrections and the VMD contributions into a single photon 
structure function = J2i X( li( x > Q 2 ) aim nx the quark distributions at some 
input scale Qq in analogy to the proton case. Then the photonic parton densities 
at different values of Q 2 follow from the inhomogencous evolution equations, that 
arc in LO 



l 

dq] box , a s f dz 



dt 



= h° 



(82) 



i 

dg 1 a s f dz 

~~dt ~ 2tt , 



/^[M!k + M!M ■ < 83 » 



Here, t = \n(Q 2 /A 2 ) and the inhomogeneity is given by 

box _ n 2 a r 2 , (t „,\2 



h box = 3ei-[x 2 + (l-xf] . (84) 

The solution of the homogeneous equations is similar to the solution of the 
DGLAP equations for hadrons and can therefore be called F^ had . The particular 
solution of the inhomogeneous equation is due to the inhomogeneity that stems 
from the point-like coupling of the photons to the quarks and can therefore be 
called _F 2 r ' pi . The general solution of the inhomogeneous evolution equations for 
the photon is thus given by 

F 7 = F ^had + jq.pl (85) 

This equation allows one to speak about the hadronic and the point-like part of 
the photon structure function. 

The measurement of the photon structure function is not as easy as for 
the proton case. Due to limited detector acceptance, the measured hadronic 
energy is not equal to the total hadronic energy, so that the photon energy is 
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not determined well, which leads to large systematic errors. In addition, the 
structure function is small and the cross section (78) is suppressed by 1/Q 4 , 
which leads to large statistical errors. In spite of these difficulties, the photon 
structure has been measured for various values of Q 2 . In Fig. 9 all existing data 
for F2 is presented as a function of x for different values of Q 2 (taken from [58]). 

One important difference in the behaviour of the proton and the photon 
structure functions is that F^ manifests strong scaling violation even in LO 
without gluon radiation included. It is positive in the whole x region. Further- 
more, F2 should be large at large x due to the point-like part of the photon 
structure function, while the structure function of the proton is small at large x. 
Predictions from several parametrizations for the PDF of the real photon (solid 
curve [33], dashed curve [34], dotted curve [37]) are also shown in Fig. 9. 




Fig. 9. Compilation [10] of all existing data on K7 in comparison to predictions 
of the PDF parametrizations in [33, 34, 37]. 
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4.3 Parton Distribution Functions of the Virtual Photon 

We now turn to a discussion of the virtual photon structure. Some old data 
from the PLUTO collaboration [25] exist, which show the structure function of 
a target photon with virtuality P 2 ~ 0.4 GeV 2 at Q 2 ~ 5 GeV 2 . At HERA 
information about the structure of the virtual photon has been obtained by two 
different methods [23, 24]. One was by tagging photons with a mean virtuality 
of P 2 ~ 10~ 5 GeV 2 with the electron calorimeter of the luminosity system. 
Another method was to use the beam-pipe calorimeter to tagg photons with a 
range in virtuality of 0.1 < P 2 < 0.6 GeV 2 . The ratio of the resolved to the 
direct contribution can be plotted as a function of the photon virtuality P 2 , 
where an experimental definition of the direct and the resolved part of the cross 
section has to be given. We will come back to this data in section 5, where we 
calculate one- and two-jet cross sections in NLO for ep scattering under HERA 
conditions. We will compare the ratios as defined in the experiment with our 
theoretical predictions. 

Here, we concentrate on the construction of the PDF's for a virtual photon. 
First, we state the LO QPM result for virtual photons, which substitutes the 
result (80) for real photons. The virtuality P 2 serves as a regulator in this case 
and in the limit A 2 < P 2 <C Q 2 one obtains [59, 60, 61] 

q f( x ,Q 2 )^3el^-S [ [x 2 + (l-x) 2 }\n^ + 6x(l-x)-2^ . (86) 

The evolution equations in Q 2 for virtual photons and the resulting PDF's are 
exactly calculable in perturbative QCD in a limited range of photon virtuality, 
A 2 -C P 2 <C Q 2 [60, 61]. The PDF's of the real photon are known for the region 
P 2 <C ^l 2 , as described in section 4.2. At HERA, though, the intermediate region 
p2 ^ yj2 j g Q f S p ec i a i interest, as has been noted above. The aim of Gliick, Reya 
and Stratmann (GRS) in [36] and Schuler and Sjostrand (SaS) in [37] was to 
construct PDF's for virtual photons, that are valid in the whole P 2 -region, i.e. 
< P 2 < Q 2 . We explain the constructions of these PDF's in the following. 



4.3.1 The PDF's of GRS Gliick, Reya and Stratmann have used for their 
construction a VMD inspired interpolation between the PDF's of real photons 
and those valid at P 2 > A 2 [36], since the PDF's f?(x, Q 2 ,P 2 ) obey evolution 
equations similar to those of the real photon. The question therefore reduces to 
finding appropriate boundary conditions at Q 2 = P 2 . Defining 

*^-&rfa? • (87) 

where m 2 = (0.77) 2 GeV 2 refers to some effective mass in the vector-meson 
propagator, PDF's valid for all values of P 2 are defined as 

P ' (s, q 2 = p 2 , p 2 ) = V (P 2 ).r h ' ad ( X , p 2 ) + [i - rtP 2 )]/;; (x, p 2 ) . (88) 

In this formula, the perturbatively calculable pointlike part is given by the 
function qf(x,Q 2 = P 2 ) from equation (86) and g~< *(x,Q 2 = P 2 ) = in NLO. 



Inclusive Single- and Dijet Rates in NLO QCD 



29 



For a LO construction, (x,Q 2 = P 2 ) = 0. These results are for the DIS, 
scheme [33], which is connected to the MS-scheme via the transformations 



where q]{x,Q 2 ) is given by (80). These transformations are valid only for the 
NLO distributions, whereas the LO distributions are equal in both schemes. The 
reason for introducing the DIS 7 scheme was that the differences between the LO 
and the NLO result are small. The hadronic, non-perturbative input is given by 



where p? LO = 0.25 GeV 2 and Unlo = °- 3 GcV2 - The function n{^a/ f 2 )/ 71 {x, fi 2 ) 
is just the prescription for the boundary conditions at input scale Q 2 = fi 2 for 
real photons. As one observes, 



Thus, the usual real photon PDF is regained for P 2 = 0, whereas the perturba- 
tively calculable part dominates for P 2 » A 2 . The number of flavors is set to 
Nf — 3. The heavy quark sector (c, 6, . . .) is supposed to be added as predicted 
by perturbation theory of fixed order with no active c and b quarks in the proton 
and photon PDF's. In LO this amounts to adding the processes j*g — > cc and 
7*<7 — ► bb to the cross section, keeping m c , mt, 0. 

In [27] GRS have provided PDF's of the virtual photon in a parametrized 
form in LO that can be conveniently used for numerical calculations. The input 
scale is Qo = 0.5 GcV and the restriction P 2 < Q 2 /5 is implemented as to 
fulfill the condition P 2 <C Q 2 . We show the x-distribution for the up-quark at 
a scale of Q 2 — 50 GeV 2 for three different values of P 2 , namely P 2 = 0, 1 and 
5 GeV 2 in Fig. 10 a. As one observes, the distribution decreases with increasing 
P 2 . For P 2 = the real photon PDF of Gliick, Reya, Vogt [33] is reproduced 
exactly and the curves fall on top of each other. The use of the MS scheme for 
the LO distributions has to be explained. The authors GRS and SaS both give 
distributions in the DIS 7 scheme, but their schemes differ slightly. SaS actually 
also give MS distributions. In order to make the results comparable, we treat the 
distributions of GRS formally as in NLO and use the transformation equations 
(89) also for the GRS LO distributions. 

It should be mentioned, that GRS have calculated NLO distributions in [27]. 
Distinct differences occur for larger P 2 and x > 10~ 3 which is mainly due to 
the different NLO perturbative boundary condition at P 2 = Q 2 , which does not 
exist for the real photon structure function. 



Idis,, = + 9i{x, Q 2 = mf) 

9DIS-, = 9jrs ' 



(89) 




(90) 



7j(P 2 = 0) = 1 and 1 - 7j(P 2 > yl 2 



(91) 



4.3.2 The PDF's of SaS Schulcr and Sjostrand represent the solution of the 
inhomogeneous evolution equations of the real photon as a sum of a perturbative 



Inclusive Single- and Dijet Rates in NLO QCD 



30 




Fig. 10. (a) The GRS LO prediction in the MS scheme for the up-distribution of 
a virtual photon at Q 2 — 50 GeV 2 and various fixed values of P 2 . For comparison, 
the LO prediction of GRV for a real photon [33] is shown, which lies exactly on 
top of the PDF of GRS for P 2 = 0; (b) SaSIM PDF's. 



and a non-perturbative term [37] 

Q 2 

p {x ,Q 2 ) = J2^P' VMD (x,Q 2 ;Q 2 ) + ^J22el f d *r^ M 2 -k 2 ) . 

(92) 

Here, Qq > A 2 is the input scale for the non-perturbative solution fi> VMD of the 
homogeneous evolution equations, which can be interpreted as a fluctuation of 
the real photon into vector-mesons. The second term represents the anomalous 
perturbative solutions of the 7 — > qq fluctuations, where k 2 is the virtuality of 
the gg-pair, which has a continuous spectrum. As noted above, the evolution 
equations of the PDF's of the virtual photon can be exactly calculated in the 
range Ql < P 2 <C Q 2 . For real photons in the region of P 2 <C Q%, the PDF's 
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Fig. 11. (a) Comparison between the GRS and SaSIM LO predictions for the 
up-distribution of a virtual photon at Q 2 = 50 GeV 2 and P 2 = 1.0 GeV 2 and 
the purely perturbative contribution in the MS scheme: (b) P 2 = 10 GeV 2 . The 
distribution for GRS lies exactly atop of the SaSIM curve. 



are given by equation (92). To obtain results valid for the whole P 2 region, SaS 
make use of the dispersion relation 

Q 2 

/ 7 '(^^ 2 ) = S 2 4/f(^) 2 pM^ 2 ) • (93) 

o 

This model provides the correct behaviour for both P 2 — > and the above 
described perturbative region. Now, the region of low k 2 can be associated with 
the discrete set of vector-mesons, so that by introducing the cut-off Ql in the 
/c 2 -integration SaS obtain [37] 
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+ SE^./|(i4)V-- M V) .(94) 

Qo 

These parton distributions are the solutions of the inhomogeneous evolution 
equations of the virtual photon. Note, that the input scale for the VMD PDF's 
has been shifted from Qq — > Qq with Qq < Qq. This is motivated by a study 
of the evolution equations in [63], which shows that the evolution for virtual 
photons starts later in Q 2 . 

The authors SaS have provided the PDF's of the virtual photon in a para- 
metrized form in LO for Nf = 4 in [37] for two different input scales, namely 
Qo = 0.6 GeV and Qo = 2 GeV. We will use the lower scale in this work. 
In contrast to the GRS parametrization, the c quark is included as a massless 
flavor in the PDF that undergoes the usual evolution as the other massless quarks 
except for a shift of the starting scale Qq. We show the up-quark distribution in 
comparison to the ones obtained by GRS in Fig. 10 b for the same scale and P 2 
values as in Fig. 10 a. They show roughly the same behaviour and deviate only 
in the small x region. For P 2 — 0, the GRV [33] distribution of the real photon 
is recovered more or less. We have used the SaSIM parametrization, which is 
given in the MS scheme, to make the SaS results comparable with the GRS 
distributions. 

In section 3.5 we have calculated the part of the hard cross section for an 
incoming virtual photon that couples directly to the subprocess, leading to the 
logarithmic terms ln(P 2 /Q 2 ). The logarithm is absorbed into the PDF of the 
virtual photon. As we suspect from the above discussion, the PDF of the virtual 
photon should reduce approximately to the perturbativcly calculable contribu- 
tion (72) in the region of large P 2 . Thus, comparing the contribution (72) directly 
with the, say, up-quark distribution for virtual photons should lead to results of 
the same order of magnitude for large P 2 . They cannot give exactly the same 
results, since the perturbative part of the u-quark distribution is evolved with 
help of the evolution equations. In Fig. 11 a, b we show the purely perturba- 
tive contribution in comparison to the w-quark distributions GRS and SaSIM 
at Q 2 = 50 GeV 2 for the two values P 2 = 1 GeV 2 and P 2 = 10 GeV 2 in the 
MS-scheme. As one observes, the perturbative solution and the u-distribution 
coincide rather well for the larger P 2 value, especially in the large x range, with 
a slight enhancement of the perturbative curve near x = 1 . 

5 Electron-Proton Scattering at HERA 

We come to an analysis of inclusive jet-rates in electron-proton scattering for 
slightly off-shell photons. After the introduction of the hadronic cross section, 
using the calculated partonic cross sections of section 3, we explain the matching 
of theoretical and experimental jet definitions. Afterwards, some numerical tests 
are discussed, and finally numerical results for one- and two-jet inclusive cross 
sections are given. A comparison with present HERA data is shown. 
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5.1 Hadronic Cross Section 

We write electron-proton scattering for the production of two jets as 

e(k)+P(p)^e(k')+3et 1 (E Tl , m ) + Jet 2 (E T2 , m )+X . (95) 

Here, k and p are the momenta of the incoming electron and proton, respectively, 
and k' is the momentum of the outgoing electron. The jets in the final state 
are characterized by their transverse momenta and rapidities which are 
observables in an experimental setup. The interaction of the electron with the 
proton is mediated by an electroweak vector boson with four-momentum q = 
(k — k') and virtuality P 2 = (— <? 2 ). The process is dominated by a photon, 
especially for the small virtualities under consideration. We therefore concentrate 
on the photon and neglect contributions from the other electroweak bosons. The 
phase space of the electron can be parametrized by the variables y = {pq)/{pk) 
and P 2 . In the case of small virtualities P 2 <C q 2 , where qo is the energy of the 
virtual photon, y gives the momentum fraction of the electron energy E e , carried 
away by the virtual photon, so y ~ qo/E e . The cm. energy of the hadronic 
system is given by sr = (p + k) 2 , whereas the cm. energy of the photon-proton 
subsystem is W 2 = (p + q) 2 . 

We have discussed the factorization of hard and soft regions of the electron- 
proton cross section in section 2. The hadronic cross section da H may be written 
as a convolution of the hard scattering process da e / k with the PDF of the proton 
fk/p( x b), where x b is the momentum fraction of the parton from the proton: 

da H (s H ) dx b f k/P (x h )da e/k (x b s H ) . (96) 

k J 

The hard process is given by the squared matrix element \M\ 2 , which has to be 
divided by the flux factor 2snXb and multiplied by the phase space of n final 
state particles of the subprocess and the electron, dPS^" +1 - ) : 

da e/k = ^-L-\M\ 2 dPS^ . (97) 

The matrix element \A4\ 2 separates into the hadron tensor H^" and the lepton 
tensor L^ v = 4(fc M fc£, - k'^k v - g^kk'): 

A-rrrv 

\M\ 2 = Tr L^H^ . (98) 

The constant a is the electromagnetic coupling constant. The separation of the 
phase space into a part depending only on the electron dL and a part depending 
only on the final state particles of the subprocess ciPS^™^ is easily achieved by 
inserting a delta function for the intermediate virtual photon and gives 

dPS(" +1 ) = dLdPS^ with dL = H- d ± d y^l . (99 ) 

I67H 2tt P z 

Here <j> is the azimuthal angle of the outgoing electron. This degree of freedom 
can be integrated out, yielding 

1 [d4> T l + (l-y) 2 A{l-y) + l + {l-y) 2 
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with the definitions H g = —g^ v H^ v and H L = (4P 2 )/ (s H y) 2 P^P v H^. Since we 
will consider the range of small photon virtualities P 2 throughout this work, 
the contribution Hl proportional to P 2 will be neglected. We approximate the 
spectrum of the virtual photons by 

df-y/eiv) _ a l + (l-y) 2 1 



dP 2 2tt v P 2 



(101) 



which is the unintegrated Weizsacker- Williams [8] formula. For later use we de- 
fine the virtuality P 2 ^ = 0.058 GeV 2 . By inserting P 2 ^ into the unintegrated 
Weizsacker- Williams formula, we obtain the value for the Weizsacker- Williams 
formula integrated in the region P^in — — ^max — ^ GeV 2 using the mini- 
mum photon virtuality P^ in := ■ In this way, we can reproduce the P 2 ~ 
results. 

Defining the partonic cross section 

da llk = ^-F 7/fed PS(") , (102) 

the hard scattering (97) integrated over the angle cf) can be written as 

d<J e /k = d(J l/k dj l/e {y)dy . (103) 

As discussed in the previous section, a photon with moderate virtuality in- 
teracts with a proton not only as a point-like particle, but also via its hadronic 
content. The hadronic structure of the photon is described by a PDF / 7 /;(x a ), 
introducing the new variable x a that gives the momentum fraction of the parton 
from the photon. To simplify the notation, the case of a direct photon is included 
into the PDF of the photon via the delta function / 77 (ir a ) = 5(1 — x a ). Summa- 
rizing the above results, the hadronic cross section da H (sh) may be written as 
a convolution of the hard scattering da k /i with the PDF of the photon / 7 /;(x a ) 
and of the proton f k /p(xb), multiplied by the photon spectrum df 1 / e (y): 

da H = y~] / dx a dx b dy df l/e (y) f 1 /i{x a )d<j l / k f k / P (x b ) ■ (104) 
k,l J 

The hard cross section da k i now describes the interactions of the partons from 
the photon (and the photon itself) with the partons from the proton and is given 
by the trace of the hadron tensor, multiplied with the phase space of the final 
state particles, divided by the flux factor: 

da l/k = - H l/k dPS {n) . (105) 

Ax a x h ys H 

The factorization (104) is visualized in Fig. 12. 

For electron-proton scattering, the SR case, in which the virtual photon cou- 
ples directly to the partons from the proton, and the DR case, in which the 
photon serves as a source of partons, are present. As mentioned in section 2, the 
matrix elements cannot be integrated over the whole region of phase space in 
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q 

= photon remnant 

Jets (p { ) 
proton remnant 



Fig. 12. Factorization of hard and soft contributions in electron-proton scatter- 
ing. 

NLO. A cut-off has to be introduced that separates the singular from the finite 
regions. The results for the integration over the singular regions were given in 
section 3 and the factorization of singular terms has been discussed. Although 
the calculation of the NLO matrix elements is straightforward, especially using 
algebraic programs like REDUCE [64] , the results for the full matrix elements are 
too cumbersome to be stated here. For both, the SR and the DR cross section, 
we have a set of two-body contributions and a set of three-body contributions. 
Each set is completely finite, as all singularities have been canceled or absorbed 
into PDF's. Each part depends separately on the phase space slicing parameter 
y c . The analytic calculations are valid only for very small y c , since terms 0(y c ) 
have been neglected in the analytic integrations. As explained in section 3.3, the 
two separate pieces have no physical meaning. When the two-body and three- 
body contributions are superimposed to yield a suitable inclusive cross section, 
as for example the inclusive one- or two-jet cross section, the dependence on the 
cut-off y s will cancel. This has been checked explicitly and will be demonstrated 
in section 5.4. 

We now come to the kinematics of electron-proton scattering. We first con- 
centrate on the SR cross section, for which the kinematics is most easily treated 
in the cm. system of the virtual photon and the proton, where for the three- 
vectors p + q = 0. We denote the momenta of the final state particles as p\ and 
P2, which can be expressed by their transverse momenta Et x — Et 2 — Et and 
their rapidities 771 and 772 in the 7*P cm. system by pi = ^(coshr^, 0, 0, sinhr/i) 
(remember that the azimuthal angle has been integrated out). From energy and 
momentum conservation one obtains 

W = E T {e- m +e-^) , (106) 
W 2 + P 2 

V = , (107) 

SH 

2W 

x b = 1 + £ T (sinhr7i + sinhr7 2 ) . (108) 
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The phase space, including the integration over Xb and y, can be expressed as 

dPS ^x b dy= w ^^^d m d V2 . (109) 

The Mandclstam variables s, t and u arc defined as 

s = (Pb + q) 2 = (pi +P2) 2 , 

t = (q - Pi) 2 = ( Pb - P2 ) 2 , (HO) 

U = (q - p 2 ) 2 = (p b - Pi) 2 ■ 

In the DR case, the rapidities of the final state partons and i]' 2 are expressed 
in the cm. system of the two partons and have to be boosted into the photonic 
cm. system via 

Vi = Vi + 2 ln x a for i = l,2 • (HI) 

Inserting these transformed rapidities into the above equations (106)-(109) the 
correct formulae in the case of the resolved photon for x a ,Xb,y and W, now 
containing 77-, are obtained. For x a = 1, which defines the SR case, the two 
systems are identical and 77- =77,. 

For a comparison with HERA data the rapidities and transverse momenta 
have to be transformed from the photonic cm. system to the HERA laboratory 
system. The calculation of the SR and DR cross sections proceed as for real 
photoproduction, i.e. the transverse momentum (qr) of the virtual photon and 
other small terms proportional to P 2 are neglected so that the virtual photon 
momentum is in the direction of the incoming electron and qo = E e y. The 
transformation from the cm. system into the HERA laboratory system is as for 
real photoproduction: 

,«= % + i.„^ . (112) 



5.2 Snowmass Jet Definition 

The factorization of hard and soft regions in the hadronic cross section has been 
discussed so far for the initial state. The non-perturbative and not calculable 
regions are parametrized through the PDF's of the hadron or the resolved pho- 
ton. A similar problem occurs in the final state. The partons that are emitted 
from the subprocess cannot be observed directly due to the confinement of color 
charge. The hadronization of partons into single hadrons in the final state can be 
described, similarly to hadrons in the initial state, by fragmentation functions. 
Another possibility is the observation of a shower built from a large number of 
hadrons without resolving the specific type of hadrons emitted. One then has 
to define jets in order to identify the hadron showers with individual partons or 
their combinations from the subprocess. The combination of hadrons into jets 
is done by cluster algorithms, where jet definitions can be implemented. The jet 
definitions should fulfill a number of criteria [65] , such as they should be simple 
to implement in theory and experiment, be well defined and yield finite cross 
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sections in any order of perturbation theory and give cross sections that are more 
or less insensitive to the hadronization processes. 

Looking at the theoretical side of inclusive two-jet cross sections, in LO there 
is a one-to-one correspondence between the parton from the subprocess and the 
jet in the final state. Therefore the theory is not sensitive to any specific cluster 
algorithm used in the experiment. This is not sensible, because the experimental 
results depend strongly on the used algorithm. Only in NLO can one imple- 
ment certain jet definitions on the theoretical side, because the jet can obtain a 
substructure due to the nearly collinear radiation of a parton in the final state. 

Several jet definitions have been proposed to date, one of the first being the 
(e, S) criterion of Sterman and Weinberg [66] (see also [50]). We will adopt the jet 
definition of the Snowmass meeting [67] . According to this definition, two partons 
i and j are recombined, if Ri j < R, where Ri,j — -J (r^ — r)j) 2 + {<f>i — (ftj) 2 and 
rjj, are the rapidity and the azimuthal angle of the combined jet respectively, 
defined as 



E Tj = E Tl +E T2 , (113) 
E Tl Vi + E T2 rj 2 
E Tj 

E Tl 4>i + E T2 (j> 2 



E Tx r\ x +E T2 T] 2 

VJ = 5 , (H4) 

E T , 



Etj 



(115) 



The cone-radius R is chosen as in the experimental analysis. Thus, two partons 
are considered as two separate jets or as a single jet depending on whether they 
lie outside or inside the cone with radius R around the jet momentum. In NLO, 
the final state may consist of two or three jets. The three-jet sample contains all 
three-body contributions, which do not fulfill the cone condition. The rapidities 
used for the cone constraint are evaluated in the HERA laboratory system. 



5.3 Numerical Input 

We now describe the input for the numerical calculations. We have chosen the 
CTEQ3M proton structure function [68] which is a NLO parametrization in the 
MS scheme, with vT— L = 239 MeV. This A value is also used to calculate a s 

A4 S 

from the two-loop formula 

sW {33-2N f )\n£\ (33-2^)2 m £ J y > 

We use this formula for both the LO and NLO calculations. In the case of 
photoproduction, the scale /x is set equal to the transverse momentum of the 
jets, since this is the only hard scale present in the interactions. Here, P 2 <C Ej,, 
so that we also set \i = Et- Equivalcntly, the factorization scales are chosen to 
be M 7 = M p = E T . 

For the PDF's of the virtual photon we choose either the GRS [36] set or the 
SaSIM set [37]. Both sets are given in parametrized form for all scales M 2 so that 
they can be applied without repeating the computation of the evolution. As men- 
tioned in section 4, both sets are given only in LO, i.e. the boundary conditions 
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for P 2 = M 2 and the evolution equations are in LO. Since neither of the two 
PDF's is constrained by empirical data from scattering on a virtual photon tar- 
get we consider these LO distribution functions as sufficient for our exploratory 
studies on jet production and treat them as if they were obtained in NLO. As 
noted in section 4, the heavy quarks are supposed to be added as predicted by 
fixed order perturbation theory with no active heavy quarks in the PDF's of 
the proton and the photon. Since in this section we are primarily interested in 
studying the sum of the direct and resolved contributions and the influence of 
the consistent subtractions of the NLO direct part we refrain from adding the 
LO or NLO cross sections for direct heavy quark production as suggested in 
[27, 36]. So, our investigations in connection with the GRS parametrization of 
the virtual photon PDF are for a model with three flavors only. For consistency 
we take also Nf = 3 in the NLO corrections and in the two-loop formula for 
a s . Of course, the proton PDF has been obtained for Nf = 4. In comparison to 
the GRS parametrization, we studied the relevant cross sections also with the 
virtual photon PDF's of SaS [37], which are for Nf = 4. 

The cross sections we have computed are for kinematical conditions as in 
the HERA experiments, for which positrons of E e = 27.5 GeV which collide 
with protons of E p = 820 GeV. To have the equivalent conditions as in the 
ZEUS analysis we choose the constraints y m in = 0.2 and y max = 0.8 for the 
variable occuring in the unintegrated Weizacker- Williams approximation. The 
cone radius is set to R = 1. 

5.4 Numerical Tests 

Since the separation of the two-body and three-body contributions with the 
slicing parameter y c is a purely technical device in order to distinguish the phase 
space regions where the integrations are done analytically from those where they 
are done numerically, the sum of the two- and three-body contributions should 
be independent from y c . The dependence of the two-body contributions on the 
slicing parameter is logarithmic, giving rise to (lny c )- and (In 2 y c )-terms. The 
parameter y c has to be quite small to guarantee that the approximations in the 
analytical calculations are valid. Typically, y c is of the order of 10~ 3 , forcing 
the two-body contributions to become negative, whereas the three-body cross 
sections are large and positive. In Fig. 13 a, b we have checked for two different 
values of P 2 , by varying y c between 10~ 4 and 10~ 2 , that the superimposed two- 
and three-body contributions arc independent of y c for the inclusive single-jet 
cross sections integrated over the whole kinematically allowed rj region for fixed 
Et = 20 GeV. Only the SR contribution is tested, since the insensitivity of the 
DR contributions on y c has been checked in [13, 22]. 

Furthermore, we have explicitly checked that the SR one- and two-jet cross 
sections for virtual photons arc in perfect agreement with the ones from real 
photoproduction given in [14, 13] by integrating the virtuality numerically over 
the region of small P 2 with P^ in < P 2 < 4 GeV 2 . The main contribution to the 
cross section comes from the lower integration boundary, where the dependence 
of the matrix elements on P 2 is small. 

Both, the y c -dependence test and the comparison with the results from [14, 
13], give us confidence that our computer program for the calculation of jet 
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Fig. 13. (a) Single-jet inclusive cross section integrated over the physical r\ region 
for Et = 20 GeV and for the virtuality P 2 = 0.058 GeV 2 as a function of the 
slicing parameter y c . The solid line gives the sum of the two- and the three-body 
contributions; (b) P 2 = 10.0 GeV 2 . 



cross sections in electron-proton scattering yields reliable results. It is interesting 
now to study the scale dependences of the LO and NLO cross sections. The 
relevant scales are the renormalization scale fx, the factorization scale for the 
virtual photon M 7 and that for the proton M p . Since the dependence on these 
scales should vanish in an all-order calculation, we expect the dependences to 
be reduced by going from LO to NLO. 

Of special interest is the dependence of the cross section on the factorization 
scale M 7 , which comes from the factorization of the photon initial state singu- 
larities. The dependence is logarithmic, since terms proportional to ln(M 2 /P 2 ) 
have been subtracted from the NLO cross section for the direct virtual photon, 
as indicated in equation (72). The dependence of the NLO direct single-jet in- 
clusive cross section, integrated over the region r\ G [—1.875,1.125] for Et = 7 
GeV, on the parameter M 1 /Et for two different values of P 2 is shown in Fig. 14 
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a, b as the dashed curve. It is compared to the resolved virtual photon contri- 
bution in LO, which gives the dotted curve. It is sufficient to use the LO matrix 
elements, since the main M 1 dependence of the resolved contribution stems from 
the dependence of the photon PDF on the negative logarithm — ln(M 2 /P 2 ). For 
the comparison we used the SaSIM parametrization of the virtual photon, which 
is given in the MS scheme. As one can see from the Fig. 14 a, b, the dependences 
on the logarithms of the direct and the resolved contributions cancel rather well 
in the sum. 
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Fig. 14. (a) Single-jet inclusive cross section integrated over r\ = [—1.875, 1.125] 
for Et = 7 GeV and for the virtuality P 2 = 1.0 GeV 2 as a function of the scale 
parameter M 1 /Et- The MS-SaSIM parametrization with Nj — 4 is chosen. The 
solid line gives the sum of the NLO direct and LO resolved virtual photon cross 
sections; (b) P 2 = 10.0 GeV 2 . 



Finally, we study the renormalization scale dependence of the NLO correc- 
tions as compared to the LO cross section. We consider only the SR contribu- 
tions, since the DR have been tested and shown to have a reduced scale de- 
pendence in [13, 22]. In Fig. 15 a, b we have plotted the one-jet inclusive cross 
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section integrated over r\ e [-1.875, 1.125] for E T = 7 GeV and P 2 = 1, 10 GeV 2 . 
We have used the two-loop formula for a s also in the LO calculation for better 
comparison. For P 2 = 1 GeV 2 the scale dependence of the NLO curve is reduced 
considerably. At P 2 = 10 GeV 2 though, the NLO cross section falls off slightly 
for the smaller scales below /j ~ Et- This could be attributed to the fact, that 
the scale and the virtuality are of the same order in this region and the condition 
P 2 -C Q 2 required in the construction of the virtual photon PDF begins to be- 
come violated. Actually, for a full test of the renormalization scale dependence 
it would be necessary to vary all scales \i — M 7 = M p simultaneously, since 
also the structure functions are renormalization scale dependent. It is an em- 
pirical fact, that the scale dependence of the proton structure function is small. 
The dependence of the photon structure function on the renormalization scale 
is large, but this can only be accounted for by the resolved contributions. The 
reduced renormalization scale dependence of the sum of the direct and resolved 
contributions has already been demonstrated in [13, 22]. 

5.5 Single-Jet Inclusive Cross Sections 

In this section, we present numerical results for inclusive one-jet cross sections 
as a function of the virtuality P 2 . We choose the following notation of the curves 
as to make the discussion clearer: the SR cross sections shall be denoted as Dir 
(reminding of the direct character of the virtual photon), whereas the DR cross 
sections are labeled Res. In addition, the sum of the SR and DR contributions is 
shown and labeled Sum. Note, that in the case of electron-proton scattering the D 
component does not exist. As has been calculated in section 3, large logarithmic 
contributions occur for small photon virtualities for the, direct virtual photon 
that can be subtracted and absorbed into the PDF of the virtual photon. The SR 
cross sections, where these logarithmic terms have been subtracted, arc specified 
by the index s, giving the contributions Dir s . All plots in this section are taken 
from ref. [26]. 

We first concentrate on predictions with the PDF's of GRS. In Fig. 16 a, b, c, 
the results for d 3 a/dETdr]dP 2 are shown as a function of Et integrated over r/ in 
the interval —1.125 < r/ < 1.875, which are the boundaries employed in the ZEUS 
analysis [23]. We show results for the the three values of P 2 = 0.058,0.5 and 1 
GeV 2 . For all three P 2 the cross section is dominated by the Res component at 
small Et- Near Et = 20 GeV the Dir s contribution is of the same magnitude as 
the Res cross section. The sum of the cross sections as a function of P 2 falls off 
nearly uniformly in the considered Et range with increasing P 2 . This decrease 
is stronger for smaller Et- 

Next, we studied the r) distribution of the Dir s contribution at fixed Et — 
7 GeV and the same P 2 values as in Fig. 16. The results are shown in Fig. 
17 a, b, c, where two approximations are shown, namely the LO cross section 
and the NLO cross section from [13]. There, P 2 = everywhere, except for 
the unintegrated Weizsacker- Williams approximation, which leads to a 1/P 2 
dependence. Obviously this approximation is good for P 2 = 0.058 GeV 2 . At the 
larger P 2 however it overestimates the cross section and should not be used. 
This means that the P 2 dependence of the Dir s part, although the strongest 
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Fig. 15. (a) Single-jet inclusive cross section integrated over rj = [—1.875, 1.125] 
for Et — 7 GeV and for the virtuality P 2 — 1.0 GeV 2 as a function of the scale 
parameter /j,/Et- The MS-SaSIM parametrization with Nf — 4 is chosen. The 
solid line gives the NLO direct prediction, whereas the dashed curve shows the 
LO cross section; (b) P 2 = 10.0 GeV 2 . 
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logarithmic P 2 dependence has been subtracted, should be taken into account. 
In the sum of the Dir s and Res cross sections the difference is small as long as 
the Res part dominates. This holds for the smaller Et's. The LO prediction is 
evaluated with the same structure functions and a s value as the NLO result. It is 
smaller than the NLO result, which it approaches with increasing P 2 . Of course, 
this finding depends on the chosen value of R because the NLO cross section 
depends on R, whereas the LO curve does not, as we have already noted before. 
Estimates of the inclusive cross section with LO calculations can therefore only 
be trusted for large cone radii. 

The results shown so far are for a model with three flavors only and therefore 
should not be compared to the experimental data except when the contribution 
from the charm quark is added at least in LO. A more realistic approach is to use 
the photon PDF's SaSIM [37] which are constructed for four flavors. In Fig. 18 a, 
b, c results are presented for d 3 a / 'dExdrjclP 2 integrated over r\ e [—1.125, 1.875] 
as a function of Et for P 2 = 0.058,0.5 and 1.0 GcV 2 . We can compare these 
curves with the results in Fig. 16a, b, c obtained with the PDF of GRS, where 
Nf = 3. The sum of the Dir s and Res contributions changes by 10% to 30% in 
the small Et region and approximately 50% in the large Et region. The larger 
cross section for Nf = 4 results mainly from the Dir s contribution. The direct 
component is more important for larger Et than for smaller Et- Therefore the 
increase is stronger in the large Et region. 

Of interest are also the rapidity distributions for fixed Et- These are shown 
for E T — 7 GeV as a function of rj between — 1 < rj < 2 choosing P 2 — 0.058, 1, 5 
and 9 GeV 2 in Fig. 19 a, b, c, d. We show the subtracted Dir s cross section, the 
Res cross section and their sum. The Res component has its maximum shifted to 
positive 77's in contrast to the Dir s component, as expected. The Dir., component 
decreases quite rapidly with increasing rj. This stems from the subtraction of 
the (In P 2 /M 2 ) terms as can be seen by comparison with the unsubtracted cross 
section, denoted Dir, in Fig. 19 a, b, c, d. The sum of the resolved and subtracted 
direct cross section Dir s is more or less constant for the smaller P 2 values and 
decreases with increasing rj for P 2 = 5 and 9 GeV 2 . 

In section 4.3 we have checked by a direct comparison, that the subtraction 
term (72) approximates the PDF of the photon rather well for large enough P 2 . 
Thus, for these large virtualities we expect the unsubtracted cross section (Dir) 
to be the correct one, rather than the sum of the subtracted direct Dir s and the 
resolved contributions, at least for small rj. The larger P 2 , the closer does the 
full direct cross section Dir approach the sum Res+Dir s , as can be observed in 
Fig. 19 a, b, c, d. As we have seen in section 4.3, there still is a deviation of the 
pure perturbative contribution from the evolved PDF in the small x region. This 
corresponds to the kinematic region of large rj, which is the forward direction of 
the proton. This deviation is evident in the Fig. 19 a, b, c, d as well; at P 2 = 9 
GeV 2 the two cross sections differ at rj = 2 by approximately 30 %. Another 
difference shows up in the backward direction at rj = —1. In this region, which 
corresponds to the region in the photon PDF where the perturbative component 
dominates, no deviation is expected. The cause might be the neglection of the 
transverse momentum qr of the virtual photon in the calculation of the Dir and 
Res cross sections, which becomes especially important for larger P 2 . Actually, 
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Fig. 16. (a) Single-jet inclusive cross section integrated over r\ e [—1.125, 1.875] 
for the virtuality P 2 = 0.058 GcV 2 . The MS-GRS parametrization with N f = 3 
is chosen. The solid line gives the sum of the subtracted direct and the resolved 
term; (b) P 2 = 0.5 GeV 2 ; (c) P 2 = 1.0 GeV 2 
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Fig. 17. (a) Single-jet inclusive cross sections for Et = 7 GeV and P 2 — 0.058 
GeV 2 . The MS-GRS parametrization with Nj — 3 is chosen. Only the SR part 
with subtraction (Dir s ) is plotted. The solid line gives the LO contribution. The 
dashed curve is the full NLO cross section, whereas the dotted curve gives the 
NLO cross section, where the NLO matrix elements have no P 2 -dependence; (b) 
P 2 = 0.5 GeV 2 ; (c) P 2 = 1.0 GeV 2 . 
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Fig. 18. (a) Single-jet inclusive cross section integrated over r\ e [—1.125, 1.875] 
for the virtuality P 2 = 0.058 GcV 2 . The MS-SaSIM parametrization with N f = 
4 is chosen. The solid line gives the sum of the subtracted direct and the resolved 
term; (b) P 2 = 0.5 GeV 2 ; (c) P 2 = 1.0 GcV 2 . 
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Fig. 19. (a) Comparisons of single-jet inclusive cross sections for Et = 7 GeV 
and the virtuality P 2 = 0.058 GeV 2 . The MS-SaSIM parametrization with N f = 
4 is chosen. The solid line gives the sum of the subtracted direct and the resolved 
term. The dash dotted curve is the direct contribution without subtraction; (b) 
P 2 = 1 GeV 2 ; (c) P 2 = 5 GeV 2 ; (c) P 2 = 9 GeV 2 . 



looking at Fig. 10 b from section 4, the purely perturbative curve, which will 
occur in the unsubtracted Dir component, overestimates the up-distribution at 
x = 1, which corresponds to the backward 77 region. So, also in this region, the 
sum of the Dir s and the Res components could still be a better estimate of the 
cross section, than the Dir component alone. 

It is clear that the resolved and the direct cross sections decrease with in- 
creasing P 2 for fixed rj and Et- It is of interest to know how the ratio of Res 
to the Dir cross section behaves as a function of P 2 . This has been analyzed in 
[26] . Apart from the fact that the ratios cannot be measured directly, we found 
a strong dependence of the ratio on the scheme chosen for the photon PDF and 
very large corrections when going from LO to NLO. As one can deduce from 
these results, it is not very sensible to compare the Dir and Res contributions 
directly. Rather one has to introduce a parameter that experimentally separates 
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Dir and Res contributions. We will introduce this parameter in the following 
section. 



5.6 Dijet Inclusive Cross Sections 

In comparison to single-jet cross sections, dijet cross sections provide a much 
stronger test of QCD, since they depend on one variable more. We will now 
present inclusive dijet cross sections d A a / dErdrjid^dP 2 as a function of P 2 . 
The variable Et is defined according to [12, 13] to be the transverse momentum 
of the measured (trigger) jet, which has rapidity r/i . The second rapidity 772 
is associated with the second jet, where the two measured jets are those with 
highest Et in the three-jet sample, i.e. Et x , Et 2 > Et 3 ■ 

In principle we could predict 77 distributions similar to those in [12]. Since 
experimental data on these distributions are not expected in the near future be- 
cause of limited statistics, we refrain from showing such plots here and present 
only the Et distributions integrated over the interval —1.125 < 771 , 772 < 1.875 
following the constraints of the ZEUS analysis [23]. The results for P 2 = 0.058, 0.5 
and 1.0 GeV 2 are shown in Fig. 20 a, b, c, where the full curve is given by the 
cross section d 4 a 7 'dETdrjid^dP 2 as a function of Et integrated over 771 and 772 
in the specified interval and for 0.2 < y < 0.8 (the plots in this section are taken 
from rcf. [26]). The functional dependence on Et does not change as a function 
of P 2 , only the absolute value of the cross section decreases with increasing P 2 . 

Furthermore we show the so-called enriched direct and resolved cross section 
in Fig. 20. These two contributions are defined with a cut on the variable x° bs , 
which is given by 

_ fy ' , (H7) 



where the sum runs over the two highest Et jets. The variable x° bs gives the 
fraction of the photon energy going into the two measured jets. It is a good 
estimate of the theoretically defined variable x a , that defines the fraction of the 
photon momentum participating in the hard interaction, see Fig. 12. For x a = 1, 
the photon couples directly to the subprocess, whereas for x a < 1 some of the 
photon energy goes into the production of a remnant jet, leading to resolved 
processes. Note, that in LO x a = x° bs . For experimental considerations, one 
defines the direct enriched contribution for x° bs > 0.75, whereas the resolved 
enriched component has x° bs < 0.75. Both enriched cross sections contain con- 
tributions from the direct and the resolved part. In Fig. 20 a, b, c the sum of the 
Dir and Res curves is equal to the full cross section d 4 a / ' dETdi]\dr\idP 2 with no 
cut on x obs . The curves in Fig. 20 are for the GRS parton distributions in the 
MS scheme. As to be expected, with increasing P 2 the full cross section is more 
and more dominated by the Dir component, in particular at the larger Et- This 
means that the cross section in x° bs < 0.75 decreases stronger with P 2 than in 
the x° bs > 0.75 region. This could be studied experimentally by measuring the 
ratio of the two cross sections as a function of P 2 for fixed Et- This has not 
been done yet. 




Fig. 20. (a) Dijet inclusive cross section integrated over 771, 772 G [—1.125, 1.875] 
for the virtuality P 2 = 0.058 GeV 2 . The MS-GRS parametrization with N f = 3 
is chosen. The solid line is the sum of the direct and the resolved contribution. 
The dashed line is the direct-enriched contribution with x° bs > 0.75 and the 
dotted curve is the resolved enriched contribution with x° bs < 0.75; (b) P 2 = 0.5 
GeV 2 ; (c) P 2 = 1.0 GeV 2 . 



Inclusive Single- and Dijet Rates in NLO QCD ... 50 



5 



4.5 r 

4 f 

3.5 r 

3 
2.5 

2 z r 
1.5 

1 r 

0.5 ^ 




0.1 



0.2 



0.3 

p 2 



0.4 



0.5 



0.6 



Fig. 21. The ratio of the resolved-enriched to the direct-enriched contributions 
as calculated in Fig. 20 a, b, c, integrated over Et x , Et 2 > 4 GeV in LO (dotted) 
and NLO (full) for the SaSIM parametrization with Nf = 4 compared with 
ZEUS data. 



Instead, the ZEUS collaboration [23] has presented data on the ratio r = 
Res/Dir, where Dir and Res refer to the enriched direct and resolved cross sec- 
tions. The ZEUS data in [23] has actually been obtained by integrating the 
transverse momenta of the two-jet cross sections over the region Et x , Et 2 > 4 
GeV and the rapidities in the range —1.125 < f?i, f?2 < 1.875 for various P 2 -bins. 
With the integration cut on the transverse momenta of the two hardest jets, the 
transverse momentum of the unobserved jet can vanish, which is not IR safe in 
NLO QCD. We therefore allow the second jet to have less than 4 GeV if the 
third unobserved particle is soft (i.e. has a transverse momentum of less than 1 
GeV) [12]. Through this procedure, a y c dependence is avoided. We calculated 
the ratio r as a function of P 2 up to P 2 = 0.6 GeV 2 and compared it with the 
ZEUS [23] data in Fig. 21 in LO (dotted curve) and NLO (full curve), using the 
SaSIM photon PDF with Nf = 4 flavors. We find quite good agreement of the 
NLO prediction with the data points for P 2 > 0.25 GeV 2 . The curve deviates 
from the data for P 2 ~ 0.2 GeV 2 , though, and even more for the point P 2 ~ 0, 
which lies about 30% above the prediction. Surely the photoproduction data is 
much more precise then the other points shown in Fig. 21. For photoproduction 
it has been shown in [12], that the measured enriched resolved component is 
larger than the predicted one for a small cut in the transverse momentum. This 
has been attributed to additional contributions from multiple interactions with 
the proton remnant jet in the resolved cross section, which have not been in- 
cluded in the NLO calculations. This underlying event contribution is reduced for 
larger E™ m and for cone radii smaller than 1. These problems must be present 
in the comparison shown here for the smallest P 2 value as well. As it seems, the 
underlying event contribution is also reduced by going to higher values of P 2 . 
This could be studied more directly by measuring rapidity distributions for the 



Inclusive Single- and Dijet Rates in NLO QCD 



51 



enriched resolved 7* sample as was done for the photoproduction case [12]. 
6 Photon-Photon Collisions 

In this section we predict inclusive jet rates for the case of photon-photon colli- 
sions in kinematic regions that will become available in LEP2 experiments. We 
first introduce the notation and kinematics for the hadronic cross section and 
then discuss some numerical results. 



6.1 Jet Production Cross Section and Kinematics 

To obtain a close correspondence between ep scattering discussed in the previous 
section and 77* scattering to be discussed here, we use similar notations. Thus, 
we deviate somewhat from the notation used in section 4. We start from electron- 
positron scattering for two-jet production, which may be written as 

e+(fci) + e-(fc 2 ) -» e + (fc' 1 ) + e-(fc^) + Jct 1 (S Tl ,r 7l ) + Jct 2 (S T2 ,7 ?2 ) + X . (118) 

We assume, that the interaction of the electrons is processed via the interaction 
of one quasi-real and one virtual photon, that are radiated by the electron and 
positron, respectively. Thus, we consider the subprocess 

ll{p)+lb(q) -]T(Jct) t +X , (119) 

with p = fci — k[ and q = k 2 — k' 2 . The electron-positron cm. energy is given 
by sh — (ki + k 2 ) 2 , whereas the 77* cm. energy is given by W 2 = (p + q) 2 . 
Since both, the real and the virtual photon, acquire some hadronic substructure 
through the resolved processes, it is not quite clear which is the probing and 
which is the probed photon. We therefore will not speak about target and probing 
photon, but simply of the real and the virtual photon. We define the virtuality 
of the real photon as Q 2 = —q 2 with Q 2 ~ and that of the virtual photon 
as P 2 = — p 2 . Next, the momentum fractions of the photon in the electron and 
positron y a and yt, have to be defined, which are given by 

Va =k^ and Vh = ^t ■ (120) 

Here, E 1 is the energy of the real photon and E e is the electron energy in the 
e + e~ cm. system. 

Considering jet production, the interaction of the photons have three different 
parts [20, 21, 22]. First in the direct (D) contribution, both photons can interact 
directly, which yields the QPM box diagram in LO. The NLO QCD corrections 
consist of the radiation of one additional gluon in the final state and the virtual 
corrections. Next, the single-resolved (SR) components have to be considered, 
which result from the hadronic structure of either of one of the photons. Since the 
resolved real and the resolved virtual photon have different hadronic structures, 
described by different PDF's, we will specify the SR contribution of a resolved 
virtual photon as SR*, whereas the single resolved real photon contribution is 
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D SR SR DR 

Fig. 22. The different components contributing in 77* scattering. 



simply denoted SR. Finally, the contribution from two resolved photons is called 
the double-resolved (DR) contribution. The different components are pictured 
in Fig. 22. 

The resolved photons are considered as sources for partons, which afterwards 
interact in a subprocess. The factorization of hard and soft regions in the e+e~ 
cross section is given by 



da(e + e — > jets) 



E 

k,l 



dx a dx b f 1 '/i{x a )d(7 k / l f kh (x b ) 

df 7 */e(ya)fy/e(yb)dyadyb ■ 



(121) 



It is written as a convolution of the PDF's of the virtual and the real photons 
f-y* /l(x a ) and fk/-y{%b), respectively, with the hard partonic cross section daij 
and the spectra of the photons, that are described by the Weizsacker- Williams 

approximation. The spectrum of the real photon is integrated over the low Q 2 

2 2 

region from Q 2 min = to Q 2 max = 4 GeV 2 , giving 



f-y/eiVb) 



1 + (i - y b f 



2tt 



In 



Vb 
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max 
i2 

min 



(122) 



The function / 7 »/ e is given by equation (101). The kinematics can be described 
most easily in the cm. system of the virtual photon and the electron that radiates 
7b. We start from the D case, for which from energy- momentum conservation 
one has 



W 

y a 
yb 



E T (e- m +e- r ' 2 ) 
W 2 +P 2 



sh 



1 



2W 



W 2 + P 2 



Et (sinh?7i + sinhr^) 



(123) 
(124) 

(125) 



Both variables y a and y b are integrated out and are included in the phase space 
for convenience, which gives 



dPS^dy a dy b 



W 2 



Et cLEt 



W 2 + P 2 2s H (2tt) 2 



dr]idr]2 



(126) 
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In the SR and DR cases, the additional variables x a and Xb have to be introduced. 
They give the momentum fractions of the partons in the resolved photons. We 
neglect the transverse momentum of the incoming virtual photon qr as we have 
done in the case of electron-proton scattering, so the partons are traveling in the 
direction of the incoming photons. Thus the energies of the partons p a and pb 
are given by x a E 1 , and XbE 7 , respectively. The rapidities are boosted by 

V'i = Vi + \ ^{x a x h ) . (127) 

In the DR ^ 1 and i(, ^ 1, whereas in the SR case only one of the 

variables x a or Xb is less then 1. 

6.2 Predictions for Inclusive Jet Rates 

We now come to a presentation of numerical results for the scattering of a virtual 
on a real photon. First we note that the factorization and renormalization scale 
dependences of the D contribution is tested indirectly with the tests done in 
section 5.4, since the matrix elements for the D case are proportional to the 
abclian color class for the subprocess 75 — > qqg which is included in the SR 
case. All tests hold for each color class separately, since e.g. the cancellation 
and factorization of singularities holds for each color class separately. The M 1 
dependence of the D and the SR* contribution compensate each other, just as 
for the SR and DR components, as tested in section 5.4. 

For producing our plots we assume kinematical conditions that will be en- 
countered at LEP2, where the photons are emitted by colliding electrons and 
positrons, both having the energy of E e = 83.25 GeV. We choose the configu- 
ration, where the virtual photon travels in the positive z-direction. We focus on 
one-jet cross sections and do not present results on dijet rates, since the studies 
here have only exploratory character. For the same reason we have used only the 
MS-GRS [36] parametrization of the photon PDF, for obtaining our results and 
do not consider the SaS PDF's [37]. We have implemented the PDF of GRS for 
both, the real and the virtual photon, since the GRS parametrization goes over 
into the GRV [33] parametrization for the real photon, when choosing Q 2 = 0. 
The real photon will be integrated over Q 2 using the Weizsacker- Williams ap- 
proximation for the region described before equation (122), whereas the virtual 
photon will have fixed P 2 -values. Because of the high cm. energies encountered 
at LEP2, we have set the number of flavors to Nf — 4, adding the contributions 
from photon-gluon fusion by fixed order perturbation theory. We took the value 
ylllL = 239 MeV for the QCD scale, which is also used in the a s two-loop for- 
mula, for which /x = Et- The factorization scales are set equal, as in the case of 
clcctron-proton scattering, with M 7 = M 7 * = Et- The Snowmass jet definition 
[67] is used as explained in section 5. 

The D and SR curves presented in the following are the NLO contributions 
for the direct virtual photon, where the large logarithm has been subtracted 
and should therefore be denoted D s and SR S in accordance with the notation in 
section 5. Since we do not present curves for the unsubtracted cross sections, we 
suppress the index s to simplify the notation. 
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Fig. 23. (a) Single-jet inclusive cross section integrated over r\ e [—2, 2] for 
the virtuality P 2 = 0.058 GeV 2 . The MS-GRS parametrization with N f = 4 
is chosen. The upper full curve is the sum of the D, SR, SR* and the DR 
components; (b) P 2 = 0.5 GeV 2 ; (c) P 2 = 1.0 GeV 2 . 
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Fig. 24. (a) Single-jet inclusive cross section as a function of T] for fixed Et = 
10 GeV and virtuality P 2 = 0.058 GeV 2 . The MS-GRS parametrization with 
Nf = 4 is chosen. The upper full curve is the sum of the D, SR, SR* and the 
DR components; (b) P 2 = 1 GeV 2 ; (c) P 2 = 5 GeV 2 ; (d) P 2 = 9 GeV 2 . 



In Fig. 23 a, b, c the Et spectra 1 for the virtualities P 2 = 0.058,0.5 and 
1.0 GeV 2 for the cross section d 3 a / ' dExdrjdP 2 are shown, integrated over the 
interval — 2 < r\ < 2, which are the boundaries being presently used at LEP1.5. 
As explained in section 5, the value P^ff — 0.058 GeV 2 is chosen as to reproduce 
the P 2 ~ case. As one can see, the SR (lower full) and SR* (dash-dotted) curves 
coincide in Fig. 23 a, where the real photon is approximated by the integrated 
Weizsacker- Williams formula and the virtual photon has the value of -P 2 // ■ ^ ne 
full cross section (upper full curve) is dominated by the DR component only 
in the small E T range for small P 2 values. For P 2 = 0.5 and 1.0 GeV 2 , the 
DR and D contributions are of the same order around Et = 4 GeV, but the 
DR component falls off quickly for the higher Et's, leaving the D component 
as the dominant contribution. This is expected, as the resolved virtual photon 

I thank T. Kleinwort for the consent in using his computer program for producing the SR* 
and DR curves. 
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is important for smaller P 2 and suppressed for the larger virtualities. For the 
same reason, the SR* contribution falls below the SR curve when going to higher 
values of P 2 (remember that the resolved real photon is not suppressed by P 2 ). 
In all curves, both SR contributions do not play an important role for the full 
cross section. Of course, all contributions decrease with increasing P 2 , so that 
the full cross section falls of with increasing P 2 . 

We turn to the ^-distribution of the single-jet cross section for fixed Et = 10 
GeV between -2 < 77 < 2 for the virtualities P 2 = 0.058,1,5 and 9 GeV 2 . As 
one sees in Fig. 24 a-d, the D and DR distributions for the lowest virtuality 
P^ff are almost symmetric, because of the identical energies of the incoming 
leptons. The SR curve falls off for negative 77, whereas the SR* component is 
suppressed for positive 77. Going to higher P 2 values, the D contribution stays 
more or less symmetric and dominates the full cross section, as we have already 
seen in Fig. 23 a, b, c for the larger Et values. The components containing 
contributions from the resolved virtual photon DR and SR* fall of in the region 
of negative 77 so that they become more and more asymmetric. This is clear, since 
we have chosen the virtual photon to be incoming from the positive ^-direction 
and the resolved virtual photon is falling off for higher virtualities. To observe 
the compensation of the ln(P 2 /M 2 ) term, subtracted from the direct virtual 
photon, with the similar but negative behavior of the resolved virtual photon, 
one has to compare the DR and SR, and the D and SR* components (see Fig. 
22). The DR and SR contributions are of the same magnitude in the negative 
77 region and the DR component is dominant for the larger 77 values, where the 
resolved photon is more important. We have observed these findings already for 
ep scattering. The same holds for the D and SR* distributions in the negative 
77 region, only here the D component is far more dominant then the SR* one in 
the whole 77 region. 

7 Summary and Outlook 

We have calculated single- and dijet inclusive jet cross sections for 7*73 and 7*7 
scattering through a consistent extension of methods used in the calculation of 77 
scattering and of photoproduction in ep scattering. The partonic cross sections 
for the considered reactions were calculated in NLO QCD using the phase-space 
slicing method, where a technical cut-off y c is introduced to separate singular and 
finite regions of phase space. The spectrum of the real photon was approximated 
with the integrated, whereas the spectrum of the virtual photon was approxi- 
mated by the differential Weizsacker- Williams formula. For the resolved virtual 
photon contribution we used the parton distributions of the virtual photon in a 
LO parametrization. 

For the hard cross section, we have in particular worked out the subtraction 
of singularities that appear when integrating the phase space over the collinear 
region of the virtual photon. Contrary to real photons, the singularity appearing 
for the virtual photon is not regulated in the dimensional regularization scheme 
but by the virtuality P 2 of the photon. This leads to a logarithm depending on 
P 2 , which is absorbed into the PDF of the virtual photon. Through this proce- 
dure, the PDF becomes scheme and scale dependent. The terms remaining in the 
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subtracted cross section are constructed in such a way, that the corresponding 
real photon term is obtained in the limit P 2 — > in the MS scheme. 

We have presented several tests of the numerical program [69] for the evalu- 
ation of the cross sections. We have shown, that the dependence on the slicing 
parameter y c vanishes when the regulated singular and finite contributions are 
added. Furthermore, the factorization scale dependences of the NLO direct and 
the LO resolved contributions cancel to a large extend. The renormalization scale 
dependence was shown to be reduced in NLO compared to LO. 

The jet cross sections for j*p scattering were computed under HERA con- 
ditions using the Snowmass jet definition. We presented distributions in the 
transverse energy and rapidity of the observed jet. For very small P 2 we found 
good numerical agreement between real and virtual photoproduction. For the 
larger P 2 values, the unsubtracted direct contribution corresponding to the case 
of deep inelastic scattering approximates the sum of the subtracted direct and 
resolved contribution rather well, at least for not too large rapidities. This is in 
accordance with the result that the perturbatively calculable subtracted term 
agrees quite well with the evolved quark distributions of the virtual photon PDF 
in the larger x range. Differences between the unsubtracted direct and the sum 
of the subtracted direct and resolved components can be attributed to small dif- 
ferences in the subtraction term and the quark distribution and to effects from 
neglecting the transverse momentum of the incoming virtual photon. 

Furthermore we have calculated distributions in the transverse energy for 
inclusive dijet cross sections. For experimental considerations the variable x° bs 
has been introduced, which is used for a separation of the direct and resolved 
contributions. The resolved part was defined for x° bs < 0.75, whereas the direct 
part was given for x° bs > 0.75. In this way, both contributions contained non- 
negligible direct and resolved parts. The sum of the enriched direct and resolved 
curves of course showed to be independent of the value of x° bs . As an application 
we have calculated the ratio of the resolved to the direct enriched cross sections, 
that could be compared to ZEUS data. The ratio shows significant NLO effects 
and is in good agreement with ZEUS data for P 2 > 0.2 GeV 2 . For smaller 
virtualities the experimental data points lie above the theoretical prediction, 
which can be attributed to, e.g., multiple scattering between the photon and 
proton remnants. 

The jet cross sections for 7*7 scattering were evaluated for conditions to be 
met at LEP2. As for ep scattering we used the Snowmass jet definition. We 
showed distributions in the transverse energy and the rapidity only for single in- 
clusive cross sections with one parametrization of the virtual photon. In contrast 
to the ep scattering, for 7*7 scattering one additional subprocess is encountered, 
which is the direct interaction of the real with the virtual photon. The singulari- 
ties of the real photon were regularized in the dimensional regularization scheme, 
whereas the virtual photon singularities have been handled as described above 
by subtracting the large logarithm. The direct component was shown to be the 
most dominant one for larger Er- The resolved virtual photon contributions were 
suppressed for larger values of P 2 due to the suppression of the virtual photon 
PDF for larger virtualities. 

Future investigations on virtual photoproduction will require more data on 
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single inclusive jet production at large transverse energy. A detailed dijet analysis 
of an infrared safe cross section such as d 4 a/dETdijidri2dP 2 , where the trans- 
verse energies of the two jets are not cut at exactly the same value, will provide 
an improved insight into the structure of the virtual photon. Furthermore, choos- 
ing a fc^-cluster-like jet definition with smaller cone radii will reduce both the 
uncertainties in the jet algorithm and in the underlying event. On the theoretical 
side, one possible improvement is the correct treatment of the transverse mo- 
mentum of the incoming photon for larger P 2 including a correct transformation 
from the photonic cm. frame to the HERA or LEP laboratory systems. For a 
consistent NLO treatment, the inclusion of NLO parton densities for the photon 
is necessary. These are, however, needed in a parametrized form and should also 
be studied in correlation with deep inelastic ej* scattering data. 
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A General Definitions 

To simplify the notation, we state some definitions here that will be used through- 
out the appendix. The Mandelstam variables s, t and u are defined in the usual 
way. The scale M 2 that appears in the virtual, initial and final state corrections 
is normally set equal to the virtuality of photon P 2 , only in the photoproduction 
limit P 2 — > we set M 2 = s. The Born terms that appear throughout this work 
are given by 




(131) 



(128) 



(129) 



(130) 



(132) 
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For the initial state corrections the plus distribution function 



R+{x,z) := 




(133) 



is needed. As the integration over z in the initial state singularities runs from 
Zmin to 1, the plus distribution function is defined as 



R+[g] 



i 

dz R(x, z)g(z) — J dz R(x, z)g(l), 
o 



(134) 



for any regular function g(z). This leads to additional terms not given here 
explicitly when (134) is transformed so that both integrals are calculated in the 
range [z min , 1]. The singular terms in the initial state corrections are proportional 
to the Altarclli-Parisi splitting functions 
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The plus functions appearing here are defined, in contrary to equation (134), in 
the limits [0, 1]. In the virtual corrections the function 
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appears [31], where £2(2) is the Dilogarithm function. In the limiting case P 2 
one finds for C 2 {x) 



( JL\ = JL -Itf ( JL 



(141) 



The square of the logarithm has two different values according to the sign of x: 
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Therefore, ones obtains the following three cases for the function L(x, y) appear- 
ing in the virtual corrections of the real photon with P 2 = 0: 

L(x,y) = ^-iln 2 (^ for x>0,y>0 , (143) 
L(x,y) = ^-W(-) fOT x<0,y<0 , (144) 



2 2 \y 

L{x,y) = ^-W(- X -) for { X>0 ' V< ° n ■ (145) 
2 V v) \x<0,y>0 1 1 



B Virtual Corrections 

In this subsection we give the explicit expressions for the virtual corrections that 
arise from the interference of the LO Born processes for 7*7 — > qq, j*q — > gq 
and 7*17 — > with the corresponding one-loop amplitudes. The expressions can 
be found in [31]. The results depend on the two-body variables s,t and u: 
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C Final State Corrections 

In the following we give the real final state corrections that appear when the 
2^3 matrix elements are integrated over the singular region of phase space. 
The expressions depend on the invariant mass cut-off y F and on the two-body 
variables s, t and u. Terms of order 0(e) have been neglected. The contributions 
i<2, . . . F 5 can be found in [31]. 
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D Initial State Corrections for Massles Partons 

Here, we state the parton initial state singularities as functions of the invariants 
s, t and u, the cut-off parameter yj and the additional variable of integration 
Zb- Again, terms of 0(e) have been neglect. The contributions 1%, ■ ■ -I§ can be 
found in [31]. 



a 2 a s Qj 



^P q ^(z b ) + JV C j(l - 2z b + 2zl) ln 



+ It) 1 - Zb 

~M 2 



Zb 



VI 



+ 2z 6 (l - z 6 ) 



(1 — £)CfT' 7 (s, t, tt) 



(155) 



7 2 b = a^Q 2 C F r 7 (s,t,u)i C* F 



+ ^(i - z 6 )( 4 + - 



In 



3 



M 2 



+ (i-^)(i + m( ^ + u)1 "% 



M 2 



21n 



-N c 

2 ° 



2z b 



2 ln M 2 "^ 



/-u 
+ 2i?+( — ,z 6 



ln f 1 - 



\M 2 

U 1 — Zb 



-u ( 1 - z b 
M*\ z b 

c-/-, i/l, * 1 , 2 _u 1 , 2 _ * 

5(1 - Zb) U ln « + 2 ln V-2 ln V 



1 - Zfc V i + u Vl z b 
1 

2 



Zb 



■21n 



-i / 1 - 06 



2 ^.l n f 1 + ^i^U-^ln(l 



1 - z b 



It = aa 2 s Q 2 



t + u yiz b ) ' 1 



M 2 \ z b 
t 1 - z h 



IJ_ P 



- Z6 

( z b) + T^-^s^gC^) 111 ( - - •'/' 



t + u yiZb 
-{t + u) 1 - Zb 
V M 2 ~ 



- 2 



I - z b 



Zb 



2 



Ty(s,u,t), 



(156) 



(157) 



Inclusive Single- and Dijet Rates in NLO QCD 



63 



= (Nf-l)I b 3 , 



n = 



+ 



aa 2 s Q 2 



2 -p, 

e 



+ 2P q ^ g (z b )\n 
+ 



— (t + u)l — Zb 



M 2 



Vi + 1 



1 



2 1 

2 



C F TJs,u,t) 

Zb J J 
./ 2 1/ tu 11 



2 h 'V 71/- 



- m+ {w^ Zb )- 2R+ {w^ Zb 



2 In 



— u / 1 — Zb 



M 2 v z b 
2z b , /. u 1 



1 — 26 



In 1 



-2 In 

Zb 



t ( 1 - z h 



t + u yiz h 
In 



M 2 v z b 
2z b 



1 



In 1 



A(l + zl) l - Zh --M* + «)l- 



z b 



M 2 



z b 



(158) 



2N£\ 
3N C J 



t l-z b 
t + u yiZb 



-J- ) L "•' 



2n 2 



(159) 



E Initial State Corrections for the Virtual Photon 



A virtual photon can decay into a qq-paiv. After the integration over the collinear 
region of phase space the Born matrix elements factorize. There are three types 
of Born matrix elements, namely the processes qq' —> qq' (denoted by T\), qq — > 
qq (denoted by T\ and T 2 ) and qq — > gg (denoted by T3). The divergence is 
regularized by the virtuality of the photon P 2 . The corrections depend on the 
two-body variables s, t and u, on the cut-off parameter yj and on the additional 
variable of integration z a . 
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